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Abstract 

The  Training  Range  Environmental  Evaluation  and  Characterization  Sys¬ 
tem  (TREECS™)  was  applied  to  a  groundwater  site  and  a  surface  water 
site  to  evaluate  the  sensitivity  of  the  uncertain  input  properties  (degrada¬ 
tion  rates  and  water-to-soil  and  water-to-sediment  adsorption  partitioning 
distribution  coefficients,  or  Kd  values)  on  the  fate  of  three  emerging  con¬ 
stituents  (EC):  DNAN,  NTO,  and  NQ.  These  three  EC  are  components  of 
the  new  insensitive  munitions  explosive  formulation  IMX-101.  The  au¬ 
thors  reviewed  the  literature  and  provided  estimates  of  the  expected  deg¬ 
radation  rates  and  Kd  values  for  the  three  EC.  EC  simulations  were 
conducted  for  each  site,  and  predicted  concentrations  in  receiving 
groundwater  and  surface  water  were  studied  to  draw  conclusions  regard¬ 
ing  the  sensitivity  of  the  results  to  these  uncertain  inputs. 

The  applications  revealed  that,  in  general,  surface  water  sites  require 
much  less  accuracy  for  estimating  soil  and  surface  water  EC  degradation 
rates  than  do  vadose  zone  and  groundwater.  Given  the  expected  low  Kd 
values  of  the  three  EC,  there  is  minimal  need  to  accurately  refine  the  EC  Kd 
estimates,  particularly  for  surface  water  receptors.  The  authors  concluded 
that  future  research  should  focus  predominantly  on  determining  degrada¬ 
tion  rates  of  EC  in  groundwater,  vadose  zone,  and  surface  soil,  in  that  or¬ 
der  of  priority. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
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All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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1  Introduction 

1.1  Background 

The  Training  Range  Environmental  Evaluation  and  Characterization  Sys¬ 
tem  (TREECS™  0  was  developed  for  the  Army  to  forecast  the  fate  of  and 
risk  from  munitions  constituents  (MC),  such  as  high  explosives  (HE)  and 
metals,  within  and  transported  from  firing/training  ranges  to  surface  wa¬ 
ter  and  groundwater.  TREECS™  consists  of  time-varying  contaminant 
fate/transport  models  for  soil,  vadose  zone,  groundwater,  and  surface  wa¬ 
ter;  these  models  forecast  MC  export  from  ranges  and  the  resulting  con¬ 
centrations  in  each  medium.  TREECS™  allows  Department  of  Defense 
(DoD)  training  range  managers  and/or  their  environmental  specialists  to 
rapidly  assess  off-site  MC  migration  as  well  as  other  contaminant  migra¬ 
tion  to  determine  whether  and  when  range  operations  could  pose  risks  to 
human  and  ecological  receptors  down-gradient  of  ranges.  Additionally, 
TREECS™  can  be  used  to  evaluate  Green  Range  Best  Management  Prac¬ 
tice  (BMP)  alternatives  where  concentrations  are  presently  or  are  predict¬ 
ed  in  the  future  to  exceed  protective  action  limits  for  human  and  ecological 
health.  A  detailed  description  of  TREECS™,  as  well  as  its  performance,  is 
provided  by  Dortch  et  al.  (2013). 

The  physicochemical  properties  of  the  MC  or  contaminants  (including 
emerging  contaminants,  or  EC)  of  interest  are  required  for  TREECS™  ap¬ 
plication.  Such  properties  include,  for  example,  molecular  weight,  solubili¬ 
ty,  solid-phase  density,  sorption  partitioning  coefficients,  Henry’s  Law 
constant,  and  degradation  rates.  Although  TREECS™  contains  three  sepa¬ 
rate  databases  for  constituent  physicochemical  properties,  there  are  data 
gaps  within  these  databases. 

The  Chemical  Transformation  Simulator  (CTS),  previously  called  the  Envi¬ 
ronmental  Fate  Simulator1 2  (EFS),  was  developed  by  the  US  Environmental 
Protection  Agency  (EPA)  to  provide  physicochemical  properties  of  com- 


1http://www.erdc.usace.armv.mil/Media/FactSheets/FactSheetArticleView/tabid/9254/Article/476659 

/training-range-envircnmental-evaliiaticn-characterizaticn-svstem.aspx 


2  http://www.epa.gpv/athens/research/efs  indepth.html 
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plex  organic  chemicals.  CTS  has  capabilities  for  estimating  properties  in 
the  absence  of  experimentally  obtained  properties;  thus,  CTS  can  help  fill 
data  gaps  for  properties,  particularly  for  emerging  contaminants  with  lim¬ 
ited  experimental  data. 

The  Environmental  Security  Technology  Certification  Program  (ESTCP) 
funded  the  work  presented  in  this  report.  The  work  focused  on  applying 
TREECS™  (with  supplemental  information  provided  by  CTS)  to  DoD 
study  sites  to  validate  the  capability  to  predict  MC  concentrations  in  re¬ 
ceiving  waters  down-gradient  of  training/firing  ranges.  The  effort  also 
demonstrated  the  utility  of  the  modeling  systems  for  forecasting  the  fate  of 
MC  —  as  well  as  of  EC  —  within  and  off-site  of  DoD  installations.  The  bulk 
of  this  work  is  reported  in  a  separate  report  (Dortch  et  al.  2017).  As  a  side 
benefit  of  this  work,  the  present  report  documents  the  utility  of  TREECS™ 
for  evaluating  EC  fate  and  provides  important  information  regarding  the 
sensitivity  of  EC  input  parameters  and  the  effects  of  their  uncertainty  on 
field  concentrations. 

1.2  Objectives 

The  objectives  of  the  study  reported  herein  were  to  use  TREECS™  to  eval¬ 
uate  the  sensitivity  of  input  properties  for  emerging  constituent  compo¬ 
nents  making  up  the  new  insensitive  munitions  (IM)  explosive 
formulation  IMX-101  and  to  assess  the  effects  of  property  uncertainty  on 
the  fate  of  those  EC  in  receiving  water  down-gradient  of  usage  areas. 

1.3  Approach 

Three  sites  were  selected  for  this  ESTCP  project  as  detailed  in  the  ESTCP 
project  Site  Selection  Memo.  The  three  study  sites  are:  Demolition  (Demo) 
Area  2  of  Camp  Edwards,  Massachusetts  Military  Reservation  (MMR), 

MA;  the  Artillery  Impact  Area  (AIA)  of  the  U.S.  Military  Academy 
(USMA),  West  Point,  NY;  and  the  Zulu  Impact  Area  (ZIA)  of  Marine  Corps 
Base  (MCB)  Camp  Pendleton,  CA.  Two  of  these  sites,  Demo  Area  2  of 
MMR  and  the  AIA  of  USMA  (West  Point),  were  used  for  further  study  of 
IMX-101  components  as  presented  in  this  report.  Surface  soil,  vadose 
zone,  and  groundwater  were  modeled  for  Demo  Area  2  of  MMR.  Surface 
soil  and  surface  water  (Popolopen  Brook)  were  modeled  for  the  AIA  at 
West  Point.  Thus,  these  two  sites  provide  both  types  of  receiving  water 
(groundwater  and  surface  water)  that  can  exist  down-gradient  of  firing 
ranges. 
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Insensitive  munitions  (IM)  contain  explosive  constituents  that  are  less 
sensitive  to  heat  and  shock.  The  IM  explosive  formulation  referred  to  as 
IMX-101  is  featured  in  the  present  study.  IMX-101  contains  three  EC:  2,  4- 
dinitroanisole  (DNAN);  3-nitro-i,  2,  4-triazol-5-one  (NTO);  and  ni- 
troguanidine  (NQ).  Thus,  these  three  EC  were  modeled  with  TREECS™  for 
the  two  study  sites,  Demo  Area  2  of  MMR,  and  AIA  of  USMA. 

As  stated  previously,  physicochemcial  properties  must  be  provided  for 
each  explosive  constituent  for  modeling  the  constituent  fate  in  various 
media.  There  is  relatively  high  confidence  in  values  for  some  of  the  proper¬ 
ties,  including  molecular  weight,  water  solubility,  and  solid  phase  density 
for  the  three  EC.  Although  there  is  some  uncertainty  in  the  values  of  Hen¬ 
ry’s  Law  Constant  (HLC),  values  for  the  three  EC  are  low  enough  that  vo¬ 
latilization  is  a  minor  fate  process.  Likewise,  diffusivity  in  air,  which  is  an 
input  used  to  compute  volatilization,  has  low  importance  for  the  three  EC. 
The  molecular  diffusivity  in  water  is  required,  but  values  can  be  readily  es¬ 
timated  (Hayduk  and  Laudie  1974). 

The  organic  carbon  normalized  soil-water  sorption  partition  coefficient, 
Koc,  L/kg,  and  the  octanol -water  sorption  partition  coefficient,  Kow, 
mL/mL,  are  used  to  estimate  the  constituent  sorption  partitioning  distri¬ 
bution  coefficients,  Kd  (L/kg),  for  soil  -  water  and  sediment  -  water,  re¬ 
spectively.  However,  Kd  values  for  some  chemicals  depend  more  on  other 
media  characteristics,  such  as  pH  for  NTO,  rather  than  organic  carbon, 
which  relates  to  the  need  for  Koc  and/or  Kow  values.  Thus,  the  true  need  is 
an  estimate  of  constituent  Kd  values  for  soil  and  sediment.  The  remaining 
property  required  for  model  input  is  usually  the  most  important  and  the 
poorest  known,  and  it  is  the  constituent  degradation  rate  or  half-life  in 
each  medium  (soil,  groundwater,  and  surface  water).  In  summary,  the  two 
types  of  inputs  that  are  quite  important  and  the  most  uncertain  are  media 
Kd  values  and  degradation  rates  (or  half-lives).  Thus,  sensitivity  testing 
was  conducted  for  these  two  inputs  for  each  EC  in  each  medium  to  evalu¬ 
ate  their  importance  with  regard  to  their  fate  and  concentration  in  receiv¬ 
ing  waters. 

All  site-specific  inputs  required  for  modeling,  such  as  soil  characteristics, 
meteorology,  and  hydrology,  were  the  same  as  those  presented  by  Dortch 
et  al.  (2017)  for  each  study  site.  The  source  loading  rate  (grams/year)  for 
each  EC  was  set  to  the  same  value  as  that  determined  for  the  explosive 
component  RDX  for  each  study  site.  The  General  Soil  Source  Zone  option 
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in  the  Operational  Inputs/Site  Conditions  screen  of  TREECS™  was  used 
for  specifying  the  site  loading  rate  for  each  EC. 

It  was  necessary  to  establish  input  values  for  all  physicochemical  proper¬ 
ties  for  the  three  EC  prior  to  conducting  any  modeling.  Thus,  the  next 
chapter  of  this  report  is  devoted  to  the  discussion  of  these  inputs. 
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2  Physicochemical  Properties 

The  discussion  of  values  for  physicochemical  properties  is  divided  into 
three  categories:  general  chemical  properties,  partition  coefficients,  and 
half-lives.  The  general  chemical  properties  are  held  constant  for  each  ap¬ 
plication  presented  herein,  but  the  partition  coefficients  and  half-lives  are 
varied  due  to  their  much  greater  uncertainty. 

2.1  General  chemical  properties 

The  general  chemical  properties  are  divided  here  into  two  groups,  those 
that  are  well-known  basic  properties,  and  those  that  are  fairly  well-known 
with  a  range  of  values.  The  well-known,  basic  properties  include  molecular 
formula,  molecular  weight,  molecular  diffusivity  in  water,  and  solid  phase 
density.  Values  of  these  properties  for  each  EC  are  shown  in  Table  l.  The 
molecular  diffusivity  in  water  was  computed  by  the  method  of  Hayduk  and 
Laudie  (1974).  The  molecular  diffusivity  in  air  was  computed  by  the  meth¬ 
od  of  Fuller  et  al.  (1966).  The  Chemical  Abstract  Service  Registry  Number 
(CASRN)  for  each  EC  is  also  shown  in  Table  1.  All  values  in  Table  1  were 
obtained  from  the  TREECS™  Army  Range  Constituent  Database 
(ARCDB),  which  provides  the  source  references  for  each  value.  Many  of 
the  values  within  ARCDB  were  obtained  from  the  EPI  Suite  software  de¬ 
veloped  by  the  US  Environmental  Protection  Agency’s  Office  of  Pollution 
Prevention  Toxics  and  Syracuse  Research  Corporation 
f http://www.epa.gov/oppt/exposure/pubs/episuitedl.htm).  EPI  Suite  is  also  part  of  CTS. 


Table  1.  Basic  properties  for  the  three  EC. 


Emerging 

Constituent 

CASRN 

Molecular 

formula 

Molecular 

weight 

Molecular 
diffusivity 
in  water, 
cm2/sec 

Molecular 
diffusivity 
in  air, 
cm2/sec 

Solid 

phase 

density, 

^cm3 

2,4- 

din  itroan  iso  1 
(DNAN) 

119-27-7 

C7H6N2O5 

198 

6.2E-6 

0.065 

1.34 

nitrotiazolone 

(NTO) 

932-64-9 

C2H2N4O3 

130 

8.5E-6 

0.1 

1.93 

nitroguanidine 

(NQ) 

556-88-7 

CH4N4O2 

104 

8.4E-6 

0.1 

1.72 
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The  fairly  well-known  chemical  properties  required  within  TREECS™  in¬ 
clude  water  solubility,  HLC,  Kow,  and  Koc,  the  latter  two  properties  are  only 
needed  for  estimating  respectively  the  soil  and  sediment  partitioning  coef¬ 
ficients,  which  are  addressed  directly  in  this  chapter.  The  ARCDB  provides 
values  for  all  four  of  these  properties,  and  in  some  cases,  multiple  values 
are  provided  so  that  the  user  must  select  which  one  to  use  in  the  applica¬ 
tion.  Source  references  are  provided  for  each  value.  The  recommended 
values  used  for  the  applications  presented  in  this  report  are  shown  in  Ta¬ 
ble  2. 


Table  2.  Fairly  well-known  chemical  properties  for  the  three  EC. 


Emerging 

Constituent 

Solubility, 

mg/L 

HLC, 

atm- 

m3/mole 

Kow , 

mL/mL 

Koc, 

L/kg 

2,4- 

dinitroanisol 

(DNAN) 

276 

3.01E-7 

43.6 

230 

nitrotiazolone 

(NTO) 

16,600 

4.07E-13 

2.34 

0.41 

nitroguanidine 

(NQ) 

2,600 

4.49E-12 

0.13 

12 

2.2  Partition  coefficients 

There  are  two  classes  of  partition  coefficients  required  for  the  two  applica¬ 
tions  presented  herein,  soil-water  and  sediment -water.  The  reversible, 
equilibrium,  linear,  partitioning  distribution  coefficient,  Kd  (L/kg),  is  used 
in  the  TREECS™  fate  models  for  various  media.  Values  for  EC  partition 
coefficients  for  the  two  classes,  soil  and  sediment,  are  discussed  in  the  re¬ 
spective  sections  below. 

2.2.1  Partition  coefficients  for  soil 

Soil  can  be  further  divided  into  surficial  soil  and  deeper  soil,  including  the 
vadose  zone  and  groundwater.  Thus,  different  Kd  values  can  be  required 
for  each  of  these  three  medium,  but  the  variation  is  usually  due  to  a  char¬ 
acteristic  of  the  medium,  such  as  the  organic  matter  content,  the  pH,  or 
some  other  medium  feature.  Sediments  refer  to  solids  within  the  benthic 
bed  and  those  that  are  suspended  in  the  water  column  of  surface  waters. 
Thus,  different  Kd  values  can  be  required  for  benthic  and  suspended  sedi- 
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ments,  and  the  variation  is  usually  due  to  a  characteristic  of  the  sediments, 
such  as  the  organic  carbon  content  and  the  solids  concentration.  Partition 
coefficients  can  vary  widely  for  a  specific  chemical  and  can  be  very  site- 
specific.  Thus,  there  can  be  much  uncertainty  associated  with  their  values. 

Scientific  literature  was  used  to  obtain  estimates  of  Kd  values  in  soils  for 
the  three  EC.  No  information  was  found  for  Kd  values  in  sediments  for  the 
three  EC,  so  sediment  values  had  to  be  extended  from  soil  values. 

The  most  comprehensive  study  of  partitioning  of  DNAN  and  NTO  in  soil 
was  conducted  by  Dontsova  et  al.  (2014a).  Empirical  relationships  for  es¬ 
timating  Kd  values  for  DNAN  and  NTO  in  soils  are  provided  by  Dontsova 
et  al.  (2014b),  where  Kd  (L/kg)  for  DNAN  depends  on  the  soil  percent  of 
organic  carbon  (OC),  and  Kd  (L/kg)  for  NTO  depends  on  soil  pH.  These 
two  relationships  for  Kd  values  of  DNAN  and  NTO  are 

Kd(DNAN)  =  1.0  +  1.  IOC  (1) 

Kd(NTO)  =  1.07-0. 12  pH  (2) 

TREECS™  has  a  tool  for  estimating  soil  Kd  values  given  Koc,  the  soil  tex¬ 
ture  (percent  sand,  silt,  and  clay),  and  the  percent  organic  matter  content. 
However,  the  predicted  DNAN  Kd  values  using  this  tool  were  about  double 
the  measured  values  (average  relative  error  of  116  %)  for  the  11  soils  stud¬ 
ied  by  Dontsova  et  al.  (2014a).  Use  of  Equation  1  provides  much  greater 
accuracy  for  estimating  soil  Kd  values  for  DNAN,  with  an  average  relative 
error  of  20  %. 

The  soil  OC  data  presented  by  Dontsova  et  al.  (2014a)  for  11  different  soils 
was  used  with  Equation  1  to  compute  DNAN  soil  Kd  values.  The  measured 
Kd  values  were  subtracted  from  the  computed  values  and  divided  by  the 
computed  value  to  produce  a  relative  difference  for  each  soil.  Then  the 
maximum  and  minimum  relative  differences  ( RD )  for  the  11  soils  were  re¬ 
tained  to  estimate  the  range  in  Kd  values  for  a  given  OC.  The  maximum 
and  minimum  RD  in  DNAN  Kd  values  for  the  11  soils  were  0.56  and  -0.56, 
respectively.  The  Kd  values  for  each  of  the  two  study  sites  were  determined 
using  the  soil  OC  values  for  each  site.  The  range  in  site  DNAN  Kd  values  for 
soil  was  then  estimated  by  adding  to  the  computed  site  Kd  value  the  prod¬ 
ucts  of  the  maximum  and  minimum  RD  times  the  computed  site  Kd  value 
as  follows, 
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Kd  (max, min)  —  Kd  +  RD(max,min)  Kd  (3) 

where  the  Kd  values  in  Equation  3  are  computed  from  Equation  1  using  the 
site  OC,  and  RD(mcoc,min)  is  the  maximum  (or  minimum)  RD  for  the  11 
soils. 

The  same  approach  described  above  for  estimating  the  range  of  DNAN  soil 
Kd  values  was  used  for  NTO,  except  that  the  soil  pH  for  the  11  soils  was 
used  with  Equation  2.  The  maximum  and  minimum  relative  differences  in 
NTO  Kd  values  for  the  11  soils  were  0.76  and  -0.3,  respectively.  The  range 
in  site  NTO  Kd  values  for  soil  was  then  estimated  by  adding  to  the  comput¬ 
ed  site  Kd  value  (using  Equation  2)  the  products  of  the  maximum  and  min¬ 
imum  RD  times  the  computed  site  Kd  value  (see  Equation  3). 

With  the  measured  soil  OC  and  pH  values  for  the  two  study  sites  as  report¬ 
ed  by  Dortch  et  al.  (2017),  the  recommended  mean,  maximum,  and  mini¬ 
mum  Kd  values  for  surficial  and  deeper  soils  were  estimated  using 
Equations  1-3.  The  surface  and  subsurface  soil  OC  at  Demo  Area  2, 

MMR,  are  0.68  and  0.068  %,  respectively,  and  the  surface  and  subsurface 
pH  values  are  4.6  and  5.9,  respectively.  The  surface  soil  OC  and  pH  at  AIA, 
USMA,  are  1.8  %  and  5.5,  respectively.  Subsurface  soil  values  were  not  re¬ 
quired  since  groundwater  was  not  modeled  at  USMA.  The  resulting  Kd 
values  are  summarized  in  Table  3. 


Table  3.  Estimated  soil  Kd  values  (L/kg)  for  two  study  sites. 


Site/ EC 

Mean 
surface 
soil  Kd 

Maximum 
surface 
soil  Kd 

Minimum 
surface 
soil  Kd 

Mean 

subsurface 
soil  Kd 

Maximum 
subsurface 
soil  Kd 

Minimum 
subsurface 
soil  Kd 

Demo  Area  2,  MMR 

DNAN 

1.75 

2.73 

0.78 

1.07 

1.68 

0.48 

NTO 

0.52 

0.91 

0.36 

0.36 

0.64 

0.25 

NQ 

0.3 

0.6 

0.1 

0.3 

0.6 

0.1 

AIA,  USMA 

DNAN 

2.98 

4.66 

1.33 

NA 

NA 

NA 

NTO 

0.41 

0.72 

0.29 

NA 

NA 

NA 

NQ 

0.3 

0.6 

0.1 

NA 

NA 

NA 

NA  is  not  applicable 
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Taylor  et  al.  (2012)  found  that  soil  Kd  values  for  NQ  varied  from  0.24  to 
0.6  L/kg  for  three  soils.  These  studies  showed  no  correlation  to  OC.  Haag 
et  al.  (1990)  stated  that  Kd  values  were  less  than  0.1  L/kg  for  two  soil 
measurements  they  conducted.  In  studies  by  Brannon  et  al.  (2004),  NQ 
partition  coefficients  were  very  low,  ranging  from  0.15  to  0.43  L/Kg.  It  is 
concluded  that  Kd  values  of  NQ  in  soil  vary  from  about  0.1  to  0.6,  with  a 
recommended  mean  value  of  about  0.3  L/kg. 

2.2.2  Partition  coefficients  for  sediments 

The  partition  coefficients  for  DNAN  and  NTO  in  surface  water  sediments 
were  estimated  based  upon  the  same  approach  used  for  soils;  i.e.,  Equa¬ 
tions  1-3.  The  total  suspended  sediment  (TSS)  and  benthic  sediment  OC 
in  Popolopen  Brook  down-gradient  of  the  AIA,  USMA,  were  both  assumed 
to  be  2  %,  which  is  close  to  the  soil  OC  and  is  a  typical  value  for  surface  wa¬ 
ter  sediments.  No  measured  values  for  sediment  OC  in  Popolopen  Brook 
were  available.  The  pH  in  Popolopen  Brook  as  reported  by  ATC  (2004)  is 
7.2.  These  values  were  used  to  compute  surface  water  sediment  DNAN  and 
NTO  Kd  values  as  summarized  in  Table  4. 

Table  4.  Estimated  surface  water  TSS  and  benthic  sediment 
Kd  values  (L/kg)  for  Popolopen  Brook,  USMA. 


EC 

Mean  Kd 

Maximum 

Kd 

Minimum 

Kd 

DNAN 

3.2 

5.0 

1.42 

NTO 

0.21 

0.36 

0.14 

NQ 

0.3 

0.6 

0.1 

The  same  Kd  values  of  NQ  in  soil  were  used  for  surface  water  sediments; 
i.e.,  a  recommended  mean  value  of  0.3  L/kg  with  a  range  of  from  0.1  to  0.6 
L/kg. 

2.3  Half-lives 

As  with  partition  coefficients,  the  scientific  literature  was  used  to  obtain 
estimates  of  degradation  rates,  or  half-lives,  for  DNAN,  NTO,  and  NQ  in 
soil  and  water.  Most  of  the  available  information  was  for  soil,  however.  In¬ 
formation  obtained  for  half-lives  of  each  of  the  three  EC  are  discussed  in 
the  respective  sections  below. 
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2.3.1  DNAN 

Half-life  of  DNAN  due  to  photo-transformation,  via  mostly  photo¬ 
oxidation,  is  about  0.7  days  (Rao  et  al.  2013).  Thus,  photo-oxidation  can  be 
a  major  loss  mechanism  in  surface  water.  Alkaline  hydrolysis  can  also 
readily  degrade  DNAN  (Salter-Blanc  et  al.  2013). 

Degradation  half-lives  for  DNAN  (as  well  as  NTO)  in  soil  on  the  order  of  a 
few  days  to  a  month  or  two  were  reported  by  Dontsova  et  al.  (2014a)  for 
their  batch  laboratory  studies.  These  results  seem  to  be  in  contradiction  to 
the  slow  rate  of  degradation  of  DNAN  in  soil  presented  by  Perreault  et  al. 
(2012a)  for  an  unamended  microcosm,  where  a  half-life  of  about  0.5  years 
was  estimated.  This  longer  half-life  is  probably  more  representative  of 
degradation  in  natural  aerobic  soil  and  groundwater  environments,  such 
as  Demo  Area  2. 

Field  et  al.  (2014)  reported  the  biotransformation  of  DNAN  and  NTO  in 
soil  cultures.  Each  soil  type  was  tested  under  aerobic  incubations  without 
any  supplemental  substrate  and  under  anaerobic  conditions  with  H2  as  an 
electron  donating  substrate.  Neither  DNAN  nor  NTO  were  significantly 
degraded  in  any  soil  under  aerobic  conditions.  However,  there  was  a  slight 
decrease  in  DNAN  concentration  over  60  days,  indicating  a  half-life  of 
roughly  2.2  years.  Both  DNAN  and  NTO  were  readily  biotransformed  un¬ 
der  anaerobic  conditions.  Anaerobic  degradation  of  DNAN  occurred  with 
half-lives  between  about  2  to  4  days.  These  soils  had  rather  high  OC. 

A  study  by  the  Environmental  Analytical  Chemistry  and  Applied  Ecotoxi- 
cology  Groups  (2012)  found  that  soil  microbes  readily  transformed  DNAN 
in  aerobic  microcosms  supplemented  with  nitrogen  (NH4CI)  and  carbon 
sources  (glucose  and  succinate).  DNAN  disappeared  at  the  rate  of  1.7  ±  0.2 
nmol  h-i  g-i  and  was  completely  removed  in  8  days.  DNAN  only  slowly 
disappeared  in  unamended  microcosms  with  the  loss  of  0.58  pmols  from 
initial  concentration  of  4  pmols  after  35  days,  indicating  a  half-life  of  about 
4.4  months  or  0.37  year.  This  suggests  that  DNAN  may  persist  longer  un¬ 
der  natural  conditions. 

Perreault  et  al.  (2012a)  reported  the  aerobic  biotransformation  of  DNAN 
in  artificially  contaminated  soil  microcosms.  DNAN  was  completely  trans¬ 
formed  in  8  days  in  soil  slurries  supplemented  with  carbon  and  nitrogen 
sources.  DNAN  was  completely  transformed  in  34  days  in  slurries  supple¬ 
mented  with  carbon  alone  and  persisted  indefinitely  in  unamended  micro- 
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cosms.  However,  the  persistent  data  does  show  a  slight  decrease  in  DNAN 
over  35  days  from  4  to  3.5  micromole,  which  translates  into  a  half-life  of 
182  days,  or  about  0.5  year. 

In  conclusion,  for  aerobic,  unamended,  natural  conditions,  the  half-life  of 
DNAN  in  soil  is  expected  to  vary  in  the  range  of  a  week  to  2.2  years.  The 
best  estimate  of  expected  half-life  in  soil  is  about  0.5  years.  Half-lives  are 
expected  to  vary  with  characteristics  of  the  soil,  but  there  is  not  presently 
enough  information  to  specify  any  such  relationships.  Half-lives  of  DNAN 
in  surface  water  are  expected  to  be  very  short  (on  the  order  of  a  day)  due  to 
photooxidation,  but  depending  on  the  environmental  conditions,  such  as 
stream  shading,  water  column  light  attenuation,  and  water  depth,  half- 
lives  in  surface  water  could  be  longer. 

2.3.2  NTO 

Krzmarzick  et  al.  (2015)  were  able  to  readily  reduce  NTO  anaerobically  in 
microcosms  using  inoculated  microbial  communities  from  seven  different 
soils.  However,  NTO  was  nonbiodegradable  in  aerobic  microcosms  with  all 
seven  soil-inoculated  communities.  These  results  are  similar  to  those  for 
RDX,  where  RDX  is  degraded  under  anaerobic  conditions  but  degrades 
little  for  aerobic  conditions.  As  with  DNAN,  microbial  enrichment  tech¬ 
niques  have  been  used  to  biodegrade  NTO  rather  rapidly  for  aerobic  condi¬ 
tions.  For  example,  Richard  and  Weidhaas  (2014a)  showed  simultaneous 
aerobic  degradation  within  four  days  for  IM  components  DNAN,  NTO,  and 
NQ  in  the  explosive  formulation  IMX101  using  soil  enrichment  cultures 
involving  sludge,  soil,  and  compost. 

Data  from  Dontsova  et  al.  (2014b)  showed  some  correlation  for  NTO  decay 
rate  in  soil  versus  soil  OC,  with  rates  that  vary  from  0.001  to  0.023  per 
hour,  or  half-life  of  28.9  to  1.25  days.  Their  conclusions  stated  that  NTO 
half-lives  varied  between  1.3  and  72  days,  where  decay  rate  increases  with 
soil  OC. 

Field  et  al.  (2014)  stated  that  neither  DNAN  nor  NTO  were  significantly 
degraded  in  any  soil  under  aerobic  conditions.  However,  both  DNAN  and 
NTO  were  readily  biotransformed  under  anaerobic  conditions.  NTO  was 
readily  converted  to  3-amino-i,  2,  4-triazol-5  one  (ATO),  and  ATO  was  not 
transformed  further  in  anaerobic  conditions.  Rapid  biodegradation  of  ATO 
in  domed  soils  under  aerobic  conditions  was  confirmed  by  release  of  am¬ 
monium  and  mineralization  of  ATO  carbon  to  CO2. 
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The  Environmental  Analytical  Chemistry  and  Applied  Ecotoxicology 
Groups  (2012)  reported  that  NTO  was  removed  slowly  (25%  NTO  removed 
after  28  days)  in  aerobic  microcosms  supplemented  with  glucose  and  suc¬ 
cinate  only.  However,  NTO  persisted  in  unamended  aerobic  and  anaerobic 
microcosms. 

Dontsova  et  al.  (2014a)  found  very  rapid  phototransformation  rates  of 
NTO,  similar  to  DNAN,  with  half-lives  of  less  than  a  day.  Thus,  NTO  is  ex¬ 
pected  to  disappear  fairly  rapidly  in  surface  water. 

Based  upon  the  above  studies,  it  is  concluded  that  NTO  does  not  degrade 
or  degrades  very  slowly  in  unamended,  aerobic  soils.  It  is  difficult  to  set  a 
range  and  recommended  values  for  half-life  of  NTO  in  soil,  but  existing 
results  suggest  that  NTO  should  degrade  more  slowly  than  DNAN  and 
possibly  NQ.  For  this  study,  a  range  of  NTO  half-lives  in  soil  of  1  to  10 
years  is  used  with  a  mean,  or  recommended,  value  of  3  years.  NTO  half- 
lives  are  expected  to  vary  with  characteristics  of  the  soil,  but  there  is  not 
presently  enough  information  to  specify  any  such  relationships.  Half-lives 
of  NTO  in  surface  water  are  expected  to  be  very  short  (on  the  order  of  a 
day)  due  to  phototransformation,  but  depending  on  the  environmental 
conditions,  such  as  stream  shading,  water  column  light  attenuation,  and 
water  depth,  half-lives  in  surface  water  could  be  longer. 

2.3.3  NQ 

Haag  et  al.  (1990)  showed  that  nitroguanidine  undergoes  base-catalyzed 
hydrolysis  for  conditions  of  25  °C  and  a  pH  of  12,  but  they  also  concluded 
that  hydrolysis  at  environmental  conditions  (25  °C,  pH  less  than  9)  was 
slow,  with  half-lives  greater  than  100  days.  Their  data  indicated  that  for  a 
pH  of  7,  hydrolysis  half-life  was  1.1  years.  They  concluded  that  NQ  de¬ 
grades  very  slowly  in  unamended,  aerobic  soils,  but  anaerobic  conditions 
can  increase  NQ  degradation.  Haag  et  al.  (1990)  stated  that  the  environ¬ 
mental  fate  of  NQ  in  surface  waters  is  dominated  by  photolysis  with  sur¬ 
face  half-lives  on  the  order  of  a  day. 

NQ  in  soil  columns  was  studied  by  Kaplan  and  Kaplan  (1985)  with  contin¬ 
uous  water  flow  through  the  columns.  NQ  could  be  degraded  if  sufficient 
carbon  (glucose)  was  added  with  the  NQ  in  the  influent  water.  The  prima¬ 
ry  product  was  ammonia.  NQ  did  not  degrade  without  sufficient  carbon. 
Their  preliminary  work  also  showed  that  NQ  will  degrade  under  anaerobic 
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conditions.  Williams  et  al.  (1989)  also  concluded  that  carbon  must  be  add¬ 
ed  to  cause  degradation  of  NQ. 

Perreault  et  al.  (2012b)  reported  that  NQ  was  efficiently  degraded  in  aero¬ 
bic  soil  microcosms  when  a  carbon  source  of  glucose  and  succinate  was 
added.  However,  NQ  persisted  in  unamended  microcosms  or  under  anaer¬ 
obic  conditions.  An  aerobic  NQ-degrading  bacterium,  Variovorax  strain 
VCi,  was  isolated  from  soil  microcosms  containing  NQ  as  the  sole  nitrogen 
source.  NQ  degradation  was  inhibited  in  the  presence  of  a  more  favorable 
source  of  nitrogen. 

Richard  and  Weidhaas  (2014a)  studied  simultaneous  aerobic  degradation 
of  the  IMX-101  formulation  constituents  DNAN,  NTO,  and  NQ,  and  deg¬ 
radation  products  were  examined.  Degradation  products  over  four  days  of 
incubation  included  nitrourea;  i,2-dihydro-3H-i,2,4-triazol-3-one;  and 
2,4-dinitrophenol.  The  enrichment  culture  maximum  specific  growth  rate 
of  0.12  h-i  and  half  saturation  constant  of  288  mg  L-i  during  degradation 
of  IMX-101  as  a  sole  nitrogen  source  suggest  that  enrichment  culture 
growth  kinetics  may  closely  relate  to  those  of  other  explosive  and  nitroar- 
omatic  compounds.  All  of  DNAN  and  most  of  the  NTO  were  aerobically 
degraded  after  nine  days  when  soil  was  amended  with  carbon.  The  order  of 
degradation  was  DNAN,  NQ,  and  NTO.  These  EC  are  suspected  to  serve  as 
a  nitrogen  source  for  microbial  growth.  However,  DNAN  was  found  to 
serve  as  either  a  carbon  or  nitrogen  source.  In  another  paper  by  Richard 
and  Weidhaas  (2014b),  they  reported  results  of  degradation  of  IMX-101 
components  in  soil.  Control  plots,  where  no  plants  and  no  enrichment  cul¬ 
ture  were  added,  were  sampled  on  days  o  and  225  and  showed  a  decline  of 
83,  83,  and  95%  in  DNAN,  NTO,  and  NQ  concentrations,  respectively. 

Mulherin  et  al.  (2005)  studied  the  stability  of  NQ  in  three  moist,  unsatu¬ 
rated  soils  under  laboratory  conditions.  This  study  yielded  a  range  of  15  to 
56  days  in  the  half-life  decay  estimates  for  NQ,  depending  on  soil  organic 
carbon,  where  half-life  decreased  with  increasing  soil  organic  carbon. 

The  above-cited  studies  indicate  a  wide  variation  of  NQ  half-lives  in  soil, 
ranging  from  days  to  possibly  years.  There  is  evidence  that  long  half-lives 
can  prevail  for  unamended,  natural  soils  with  near  neutral  pH.  Given  these 
results,  a  range  of  NQ  half-life  in  soil  from  0.3  to  3  years  with  an  average 
of  1  year  was  used  for  this  study.  Half-lives  of  NTO  in  surface  water  are  ex¬ 
pected  to  be  very  short  (on  the  order  of  a  day)  due  to  photolysis,  but  de- 
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pending  on  the  environmental  conditions,  such  as  stream  shading,  water 
column  light  attenuation,  and  water  depth,  half-lives  in  surface  water 
could  be  longer  than  a  day. 

2.3.4  Summary 

The  mean  and  range  of  recommended  EC  half-lives  for  use  in  this  study 
are  summarized  in  Table  5. 


Table  5.  Recommended  half-lives  (years)  of  EC. 
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3  Demo  Area  2,  MMR,  Application 

3.1  Background  for  model  inputs 

This  chapter  presents  the  results  of  TREECS™  application  to  Demo  Area  2 
of  MMR  to  assess  the  fate  of  the  three  EC  in  groundwater  as  affected  by 
uncertain  inputs  for  soil  -  water  partitioning  and  half-life.  The  site  de¬ 
scription  and  model  inputs  for  this  study  site  were  described  previously 
(Dortch  et  al.  2007,  Dortch  2012,  Dortch  2015,  and  Dortch  et  al.  2017)  so 
are  not  repeated  here.  However,  basic  information  and  key  features  of  the 
model  inputs  regarding  this  application  are  described  below  to  provide  the 
background  required  for  discussing  model  results. 

Three  modeling  components  of  TREECS™  were  used,  including  the  Tier  2 
soil  model,  the  Multimedia  Environmental  Pollutant  Assessment  System 
(MEPAS)  vadose  zone  model,  and  the  MEPAS  aquifer  model.  The  option 
for  average  annual  hydrology,  rather  than  daily  varying  hydrology,  was 
used  in  the  soil  model  for  this  application  due  to  the  slow  progression  of 
groundwater  transport  that  can  take  decades  to  be  realized.  A  constant  EC 
residue  loading  rate  of  1,500  grams  per  year  was  applied  for  10  years.  This 
loading  rate  is  the  same  as  the  HE  RDX  loading  rate  used  in  the  previous 
applications  to  this  site,  which  was  estimated  based  on  the  types  of  demo¬ 
lition  detonation  charges  used  (Dortch  2007). 

The  model  was  run  for  100  years  using  the  adaptive  time-step  feature  and 
minimum  time  step  of  0.001  year.  Zero  initial  soil  contamination  of  EC 
was  assumed.  Surface  soil  runoff  was  assumed  to  be  zero  for  the  highly 
permeable  soils;  thus,  all  precipitation  resulted  in  evapotranspiration  and 
infiltration  into  the  soil  layer.  The  solid  and  adsorbed  phase  EC  were  as¬ 
sumed  to  be  nondegradable  (i.e.,  a  high  value  of  1E20  years  was  used  for 
half-life).  Soil  properties  for  silty  loam  and  sand  were  used  for  surface  soil 
and  vadose/aquifer,  respectively.  Chemical-specific  property  inputs  of  the 
EC  are  presented  in  Table  1  and  Table  2.  The  mean  and  range  of  recom¬ 
mended  EC  soil  Kd  values  are  presented  in  Table  3.  The  mean  and  range  of 
recommended  EC  half-lives  are  presented  in  Table  5. 

The  initial  particle  size  of  solid  phase  EC  residue  was  set  to  2791  pm  (mi¬ 
crons)  using  the  approach  recommended  by  Dortch  et  al.  (2017)  that  was 
based  on  the  modeling  of  IMX-101  fate  laboratory  studies.  This  approach 
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computes  the  mass-weighted  average  particle  size.  The  mass-weighted  av¬ 
eraged  size  was  computed  using  the  percent  by  mass  of  each  constituent  in 
IMX-101  and  the  representative  particle  size  of  each  constituent,  resulting 
in  the  input  particle  size  of  2791  pm.  This  particle  diameter  was  used  as 
the  input  value  for  each  of  the  three  EC.  The  representa-tive  particle  size  of 
6000  pm  was  used  for  DNAN  when  calculating  the  mass-weighted  size. 
DNAN  is  the  molten-cast  matrix  within  which  the  smaller  crystals  reside 
(i.e.,  NTO  and  NQ  for  IMX-101).  Thus,  the  particle  size  of  DNAN  was  as¬ 
sumed  to  be  the  same  as  that  estimated  for  RDX,  which  is  the  single  HE 
component  within  the  explosive  formulation  C4  used  at  Demo  Area  2. 

The  well  location  for  monitoring  model  output  for  groundwater  concentra¬ 
tion  was  set  to  the  location  of  monitoring  well  161  (MW161)  at  Demo  Area 
2,  which  is  one  of  the  wells  where  RDX  concentrations  were  measured  in 
the  field. 

3.2  Baseline  results 

TREECS™  was  applied  using  the  inputs  discussed  in  the  previous  section 
to  establish  baseline  results  for  each  EC.  The  recommended  mean  soil  Kd 
values  and  mean  half-lives  (HL)  in  surface  and  subsurface  soil  were  used 
for  the  baseline  application.  RDX  was  also  included  in  the  baseline  results 
for  comparing  the  EC.  The  inputs  for  RDX  were  the  same  as  those  of  the 
validation  application  presented  by  Dortch  et  al.  (2017).  The  HL  in  water- 
dissolved  phase  in  soil,  vadose,  and  aquifer  was  set  to  100  years  for  the 
modeling  of  RDX. 

The  aquifer  well  concentrations  in  parts  per  billion  (ppb),  or  micrograms 
per  liter  (pg/L),  for  the  three  EC  and  RDX  are  plotted  versus  time  as 
shown  in  Figure  1.  The  well  concentrations  of  the  three  EC  are  much  lower 
(far  below  any  reasonable  detection)  than  that  of  RDX.  The  reason  for  this 
is  due  to  the  slow  degradation  rate  (relatively  high  HL)  of  RDX,  which  was 
based  on  previous  validation  modeling  (Dortch  et  al.  2007,  Dortch  2012, 
Dortch  2015,  and  Dortch  et  al.  2017).  The  three  EC  are  almost  totally  de¬ 
graded  within  the  vadose  zone  due  to  their  much  greater  degradation  rates 
(i.e.,  smaller  HL  values)  and  the  long  transit  time  in  that  medium.  Also, 
the  three  EC  arrive  at  the  receptor  well  later  than  RDX  due  to  their  higher 
Kd  values. 
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Figure  1.  Computed  aquifer  concentration  of  EC  compared  with  RDX  down-gradient  of 
Demo  Area  2,  MMR,  for  baseline  conditions. 


Aquifer  dissolved  concentrations  at  MW161  down-gradient 
Demo  Area  2,  MMR 


of 


•RDX 
•DNAN 
•  NTO 
■NQ 


Given  the  above  baseline  results,  the  decision  was  made  to  apply  the  same 
HL  values  for  each  EC  as  used  for  RDX;  i.e.,  loo  years  for  dissolved  phase 
in  all  media,  so  that  the  fate  of  the  EC  relative  to  RDX  could  be  better 
evaluated.  The  results  of  assuming  an  HL  of  loo  years  in  soil  and  ground- 
water  for  each  EC  are  presented  in  Figure  2.  Although  the  well  peak  con¬ 
centrations  of  NTO  and  NQ  occur  later  than  that  of  RDX,  these  peak 
concentrations  are  greater  than  that  of  RDX  due  to  the  faster  dissolution 
rate  associated  with  the  much  higher  water  solubility  of  NTO  and  NQ.  Alt¬ 
hough  the  dissolution  rate  of  DNAN  is  a  little  faster  than  that  of  RDX,  its 
higher  soil  partitioning  causes  substantial  delay  and  attenuation  of  peak 
aquifer  concentration,  resulting  in  a  peak  concentration  that  is  similar  to 
that  of  RDX  for  the  same  degradation  rate,  or  HL.  Given  the  great  uncer¬ 
tainty  regarding  the  EC  degradation  rates  and  the  evidence  that  EC  degra¬ 
dation  can  be  quite  slow  in  natural,  unamended,  aerobic  groundwater,  the 
results  in  Figure  2  may  not  be  unrealistic. 

The  next  application  consisted  of  baseline  conditions  with  two  exceptions, 
soil  constituent  loadings  and  HL  values.  For  this  application,  the  constitu¬ 
ent  loading  rates  to  soil  were  held  constant  at  1500  g/yr  for  200  years,  and 
HL  values  for  all  constituents  were  set  to  1E20  years  in  all  media,  or  essen¬ 
tially  no  degradation.  The  application  was  run  for  200  years.  The  results  of 
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this  application  are  not  shown,  but  they  showed  that  RDX  and  all  three  EC 
eventually  reached  a  steady-state  aquifer  well  concentration  of  9  ppb,  but 
some  constituents  took  longer  than  others  to  reach  steady  state.  Both  NTO 
and  NQ  reached  steady  state  much  faster  than  DNAN  and  RDX  due  to 
their  high  solubility  (and  dissolution  rate)  and  low  Kd  values.  Steady  state 
was  reached  in  about  55  and  60  years  for  NQ  and  NTO,  respectively. 
DNAN  required  about  150  years  to  reach  steady  state,  while  RDX  had  not 
quite  reached  steady  state  after  200  years  due  to  its  much  lower  solubility. 


Figure  2.  Computed  aquifer  concentration  of  EC  compared  with  RDX  down-gradient  of 
Demo  Area  2,  MMR,  for  100-year  HLfor  all  constituents  in  all  media. 


3.3  Tracer  study 

A  simulation  of  a  miscible,  nonreactive,  nonadsorbing  tracer  was  conduct¬ 
ed  to  determine  the  residence  time  in  surface  soil  and  transit  time  through 
the  vadose  zone  and  aquifer  to  the  receiving  well.  This  simulation  was 
conducted  by  selecting  one  of  the  EC,  which  was  DNAN  in  this  case,  and 
setting  its  Kd  values  to  nearly  zero;  i.e.,  1E20,  and  its  HL  values  to  essen¬ 
tially  infinity  (1E20  years)  for  all  media.  The  constituent  was  set  to  misci¬ 
ble  in  the  soil  model  inputs,  which  means  that  its  solid  phase  dissolves 
instantly.  The  constituent,  or  tracer,  loading  rate  to  soil  was  set  to  a  con¬ 
stant  value  for  100  years,  and  the  simulation  was  run  for  200  years.  Tracer 
concentrations  were  plotted  versus  time  for  soil,  vadose  zone,  and  aquifer 
receptor  well  (MW161). 
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The  concentration  plot  for  soil  revealed  that  the  tracer  residence  time  in 
soil  is  approximately  three  months.  The  other  two  plots  showed  that  the 
transit  time  through  the  vadose  zone  and  aquifer  from  the  vadose  zone  en¬ 
try  to  receptor  well  are  approximately  14  years  and  2  months,  respectively. 
The  exact  transit  time  through  the  aquifer  to  the  well  was  difficult  to  de¬ 
termine  due  to  the  coarseness  of  model  output  relative  to  the  length  of  the 
simulation.  However,  with  a  Darcy  velocity  of  too  cm/day,  an  effective  po¬ 
rosity  of  0.3,  and  well  distance  of  198  m  from  the  source,  the  travel  time 
should  be  two  months  from  vadose  entry  to  receptor  well. 

The  tracer  application  revealed  that  the  transit  time  through  the  vadose 
zone  is  by  far  the  most  important  aspect  with  regard  to  constituent  degra¬ 
dation  and  fate  for  this  study  site.  However,  not  all  groundwater  sites  will 
have  this  characteristic.  For  many  sites,  transit  time  and  degradation  with¬ 
in  the  aquifer  will  be  more  important  than  for  the  vadose  zone.  This  site 
has  a  thick  vadose,  or  water-unsaturated  zone,  and  a  relatively  short  aqui¬ 
fer  travel  distance  from  vadose  entry  to  receptor  well.  For  many  sites,  the 
opposite  can  be  the  case;  i.e.,  the  vadose  zone  thickness  can  be  small  rela¬ 
tive  to  the  aquifer  transit  distance  to  the  receptor  well. 

Residence  and  transit  times  for  specific  EC  with  adsorption  are  quite  dif¬ 
ferent  from  that  associated  with  a  nonabsorbing  tracer.  Each  of  the  three 
EC  were  next  tested  individually  to  evaluate  their  residence  and  transit 
times  and  sensitivity  to  uncertain  Kd  and  HL  values. 

3.4  DNAN  sensitivity 

3.4.1  Residence  and  transit  times 

A  simulation  of  DNAN  was  conducted  to  evaluate  its  residence  time  in  sur¬ 
face  soil  and  transit  times  through  the  vadose  zone  and  aquifer  from  va¬ 
dose  entry  to  receptor  well.  This  simulation  was  conducted  in  the  same 
manner  as  described  above  for  a  general  tracer.  Thus,  DNAN  was  treated 
as  a  miscible  constituent  with  a  constant  soil  loading  rate  for  too  years 
and  no  degradation  in  any  media.  The  exception  was  that  the  Kd  values 
were  set  to  the  recommended  mean  values  for  DNAN,  or  1.75  L/kg  for  soil 
and  1.07  L/kg  for  vadose  and  aquifer  (see  Table  3). 

The  results  of  this  simulation  showed  that  the  residence  time  in  soil  is 
about  one  year,  and  the  transit  time  through  the  vadose  zone  and  from  va¬ 
dose  entry  into  the  aquifer  to  the  receiving  well  are  about  101  and  1  years, 
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respectively.  These  times  are  much  greater  than  those  associated  with  the 
nonabsorbing  tracer.  Adsorption  has  the  effect  of  extending  residence  and 
transit  times,  which  is  referred  to  as  transport  retardation.  Given  the  much 
longer  transit  time  through  the  vadose  zone,  it  is  easy  to  see  why  the  deg¬ 
radation  rate  for  the  vadose  zone  is  so  sensitive  and  important  for  this 
study  site.  The  results  of  this  application  also  reinforce  the  importance  of 
the  soil  partitioning  due  to  its  impact  on  the  effectiveness  of  degradation. 
For  example,  higher  soil  partitioning  causes  extended  transit  time,  which 
provides  more  time  for  degradation,  thus  increasing  the  importance  of  ac¬ 
curately  estimating  the  degradation  rate. 

3.4.2  Half-lives 

The  sensitivity  of  aquifer  concentrations  to  half-life  of  DNAN  in  each  me¬ 
dium  was  evaluated.  This  sensitivity  analysis  was  conducted  by  varying  the 
HL  in  one  medium  at  a  time  while  setting  the  HL  to  too  years  for  the  other 
two  media  and  then  comparing  computed  aquifer  peak  well  concentration 
to  the  peak  computed  with  HL  set  to  too  years  for  all  three  media  (see 
Figure  2).  All  Kd  values  were  set  to  baseline  values  for  these  model  runs 
(see  recommended  mean  values  in  Table  3). 

Using  the  maximum  recommended  DNAN  HL  of  2.2  years  in  soil  (see  Ta¬ 
ble  5)  and  too  years  below  ground  in  vadose  and  aquifer  only  slightly  af¬ 
fected  the  aquifer  well  concentration,  reducing  the  peak  concentration 
from  1.22  ppb  for  100-year  half-life  in  all  media  to  1.19  ppb,  or  a  2.5  %  re¬ 
duction.  Using  the  minimum  recommended  half-life  in  soil  of  0.02  years 
and  100  years  below  ground  resulted  in  reducing  well  peak  concentration 
by  almost  a  factor  4  to  0.32  ppb,  or  a  74%  reduction  from  results  with  a 
100-year  half-life  in  all  media. 

With  the  maximum  recommended  DNAN  HL  of  2.2  years  in  aquifer  (see 
Table  5)  and  100  years  in  soil  and  vadose  zone,  the  aquifer  peak  well  con¬ 
centration  was  reduced  from  1.22  ppb  for  100-year  half-life  in  all  media  to 
0.86  ppb,  or  a  30  %  reduction.  Using  the  minimum  recommended  half-life 
in  aquifer  of  0.02  years  and  100  years  for  soil  and  vadose  zone  resulted  in 
essentially  zero  concentration  in  the  aquifer.  This  result  is  due  to  the  ra¬ 
ther  long  transit  time  of  DNAN  in  the  aquifer  to  the  receptor  well,  allowing 
time  for  complete  degradation  with  the  short  HL  value. 

Using  a  half-life  of  2.2  years  in  the  vadose  zone  and  100  years  in  the  other 
two  media,  resulted  in  a  peak  well  concentration  of  3.6E-14  ppb,  or  essen- 
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tially  zero  concentration  in  the  aquifer.  The  total  loss  of  DNAN  in  the  va- 
dose  zone  is  due  to  the  much  longer  residence  time  in  the  vadose  zone  and 
the  low  HL  used  for  DNAN. 

Aquifer  concentrations  of  DNAN  are  moderately  sensitive  to  HL  in  surface 
soil,  but  are  strongly  affected  by  HL  in  vadose  and  aquifer  for  this  site. 
These  results  are  due  to  the  long  transit  times  of  DNAN  through  the  va¬ 
dose  zone  and  in  the  aquifer.  The  transit  times  are  extended  for  DNAN  rel¬ 
ative  to  a  nonadsorbing  tracer  due  to  soil  partitioning. 

3.4.3  Soil  partitioning 

The  sensitivity  of  aquifer  concentrations  to  DNAN  Kd  values  in  each  medi¬ 
um  was  evaluated.  This  sensitivity  analysis  was  conducted  by  varying  Kd  in 
one  medium  at  a  time  while  setting  the  Kd  values  to  the  recommended 
means  for  the  other  two  media  and  then  comparing  computed  aquifer 
peak  well  concentration  to  the  peak  computed  with  recommended  mean 
Kd  values  for  all  three  media  (see  Figure  2).  All  HL  values  were  set  to  100 
years  for  these  model  runs. 

Using  the  recommended  maximum  DNAN  Kd  value  in  surface  soil  of  2.73 
L/kg  (see  Table  3)  resulted  in  an  aquifer  peak  well  concentration  of  1.21 
ppb  compared  to  1.22  ppb  for  the  recommended  mean  DNAN  Kd  value  in 
surface  soil  of  1.75  L/kg.  Using  the  recommended  minimum  DNAN  Kd 
value  in  surface  soil  of  0.78  L/kg  (see  Table  3)  resulted  in  an  aquifer  peak 
well  concentration  of  1.28  ppb  compared  to  1.22  ppb  for  the  recommended 
mean  DNAN  Kd  value.  Thus,  aquifer  concentrations  of  DNAN  for  this  site 
are  fairly  insensitive  to  surface  soil  DNAN  Kd  values.  However,  it  is  recog¬ 
nized  that  the  sensitivity  to  surface  soil  Kd  values  would  increase  for  lower 
HL  in  soil. 

The  maximum  and  minimum  DNAN  Kd  values  were  applied  for  the  aquifer 
with  all  other  Kd  values  for  the  other  two  media  set  to  recommended  mean 
values  (see  Table  3).  These  simulations  showed  that  aquifer  concentrations 
of  DNAN  were  quite  close  (about  0.01  ppb  difference)  to  those  computed 
using  the  recommended  mean  value.  Thus,  aquifer  concentrations  of 
DNAN  for  this  site  are  practically  insensitive  to  aquifer  DNAN  Kd  values. 

Based  on  the  above  transit  time  analysis,  aquifer  DNAN  concentrations 
should  be  far  more  sensitive  to  the  Kd  value  for  the  vadose  zone  than  for 
soil  and  aquifer.  Using  the  maximum  recommended  value  for  DNAN  va- 
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dose  Kd  of  1.68  L/kg,  the  peak  well  concentration  for  DNAN  was  0.58  ppb 
(52%  decrease  from  1.22  ppb)  and  occurred  after  154  years  rather  than  107 
years  using  Kd  of  1.07  L/kg.  Using  the  minimum  recommended  Kd  value  of 
0.48  for  the  vadose  zone,  a  peak  well  concentration  for  DNAN  of  3.07  ppb 
(152%  increase  from  1.22  ppb)  was  computed  after  58  years.  The  percent 
of  changes  are  disproportional  due  to  the  increased  effect  of  degradation 
for  the  longer  transit  time  associated  with  the  higher  Kd  value.  The  HL  for 
all  media  was  set  to  100  years.  It  is  concluded  that  aquifer  DNAN  concen¬ 
trations  are  moderately  to  highly  sensitive  to  Kd  values  in  the  vadose  zone 
for  this  site.  Of  course,  the  sensitivity  of  DNAN  Kd  value  in  the  vadose  zone 
depends  on  DNAN  HL  in  the  vadose  zone.  Aquifer  concentration  depend¬ 
ence  on  DNAN  Kd  value  in  the  vadose  zone  is  unimportant  when  applying 
the  recommended  DNAN  HL  values  shown  in  Table  5  that  vary  between 
0.02  and  2.2  years.  Aquifer  concentrations  are  essentially  zero  for  these 
low  HL  values  regardless  of  the  DNAN  Kd  value  in  the  vadose  zone. 

3.4.4  Conclusions 

Aquifer  DNAN  concentrations  for  this  site  are  practically  insensitive  to  Kd 
values  in  soil  and  aquifer,  but  are  quite  sensitive  to  Kd  values  in  the  vadose 
zone.  The  greater  the  Kd  value  in  a  medium,  the  greater  the  effect  of  HL  in 
that  medium  on  constituent  degradation,  thus  increasing  the  need  for 
more  accurate  estimates  of  HL  for  that  medium.  Assuming  a  very  long  HL, 
aquifer  peak  concentrations  vary  with  vadose  zone  Kd  due  to  the  fact  that 
the  soil  DNAN  loading  imposed  at  this  site  extended  for  a  relatively  brief, 
finite  period  of  time  (10  years).  Had  the  loadings  been  steady  state,  steady- 
state  aquifer  concentrations  would  eventually  have  been  reached,  and  for  a 
very  long  HL,  Kd  values  would  have  had  little  effect  on  the  steady-state 
concentration  but  would  affect  the  time  to  reach  steady  state.  Aquifer  con¬ 
centrations  are  highly  dependent  on  HL  in  the  vadose  zone  due  to  the 
much  longer  transit  time  in  that  medium  for  this  study  site.  Aquifer  DNAN 
concentrations  are  dependent  on  HL  in  the  aquifer,  but  less  so  than  for 
vadose  zone.  Aquifer  concentrations  of  DNAN  are  moderately  affected  by 
HL  in  soil  if  there  is  little  or  no  degradation  in  the  vadose  zone  and  aqui¬ 
fer.  For  example,  reducing  HL  in  soil  from  100  to  0.02  years  with  a  100- 
year  HL  for  vadose  and  aquifer,  reduced  peak  well  concentration  by  a  fac¬ 
tor  of  about  4. 
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3.5  NTO  sensitivity 

3.5.1  Residence  and  transit  times 

A  simulation  of  NTO  was  conducted  to  evaluate  its  residence  time  in  sur¬ 
face  soil  and  transit  times  through  the  vadose  zone  and  aquifer  from  va- 
dose  entry  to  receptor  well.  This  simulation  was  conducted  in  the  same 
manner  as  described  above  for  a  general  tracer  and  DNAN.  Thus,  NTO  was 
treated  as  a  miscible  constituent  with  a  constant  soil  loading  rate  for  too 
years  and  no  degradation  in  any  media.  The  exception  was  that  the  Kd  val¬ 
ues  were  set  to  the  recommended  mean  values  for  NTO,  or  0.52  L/kg  for 
soil  and  0.36  L/kg  for  vadose  and  aquifer  (Table  3). 

The  results  of  this  simulation  showed  that  the  residence  time  in  soil  is 
about  10  months,  and  the  transit  time  through  the  vadose  zone  and  from 
vadose  to  the  receiving  well  are  about  41  years  and  5  months,  respectively. 
These  times  are  greater  than  those  associated  with  the  nonabsorbing  trac¬ 
er  since  adsorption  has  the  effect  of  extending  residence  and  transit  times, 
or  transport  retardation.  The  Kd  values  for  NTO  are  lower  than  those  for 
DNAN,  which  causes  faster  transit  than  DNAN  through  each  medium. 

3.5.2  Half-lives 

The  sensitivity  of  aquifer  concentrations  to  half-life  of  NTO  in  each  medi¬ 
um  was  evaluated.  This  sensitivity  analysis  was  conducted  by  varying  the 
HL  in  one  medium  at  a  time  while  setting  the  HL  to  too  years  for  the  other 
two  media  and  then  comparing  computed  aquifer  peak  well  concentration 
to  the  peak  computed  with  HL  set  to  too  years  for  all  three  media  (see 
Figure  2).  All  Kd  values  were  set  to  baseline  values  for  these  model  runs 
(see  recommended  mean  values  in  Table  3). 

Using  the  maximum  recommended  NTO  HL  of  10  years  in  soil  (see  Table 
5)  and  too  years  below  ground  in  vadose  and  aquifer  only  slightly  affected 
the  aquifer  well  concentration,  reducing  the  peak  concentration  from  4.29 
ppb  for  100-year  half-life  in  all  media  to  4.14  ppb,  or  a  3.5  %  reduction. 
Using  the  minimum  recommended  NTO  half-life  in  soil  of  1.0  year  and 
100  years  below  ground  resulted  in  reducing  well  peak  concentration  to 
4.05  ppb,  or  a  5.6  %  reduction  from  results  with  a  100-year  half-life  in  all 
media.  Thus,  aquifer  concentration  of  NTO  is  weakly  sensitive  to  NTO  HL 
in  soil  for  the  range  of  recommended  NTO  HL  values. 
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With  the  maximum  recommended  NTO  HL  of  to  years  in  aquifer  (see  Ta¬ 
ble  5)  and  too  years  in  soil  and  vadose  zone,  the  aquifer  peak  well  concen¬ 
tration  was  reduced  from  4.29  ppb  for  100-year  half-life  in  all  media  to 
4.16  ppb,  or  a  reduction  of  3  %.  Using  the  minimum  recommended  half- 
life  in  aquifer  of  1  year  and  100  years  for  soil  and  vadose  zone  resulted  in  a 
peak  well  concentration  for  NTO  of  3.08  ppb,  or  a  28  %  reduction  com¬ 
pared  with  using  an  HL  of  100  years  in  all  three  media.  Aquifer  NTO  con¬ 
centrations  are  mildly  sensitive  to  the  HL  of  NTO  in  the  aquifer  for  the 
range  of  recommended  NTO  HL  values. 

Using  an  NTO  HL  of  10  years  in  the  vadose  zone  and  100  years  in  the  other 
two  media  resulted  in  a  peak  well  concentration  of  0.39  ppb,  or  a  reduc¬ 
tion  of  91  %  compared  to  an  HL  of  100  years.  Using  a  NTO  HL  of  1  year  in 
the  vadose  zone  resulted  in  essentially  zero  aquifer  concentration  of  NTO. 
The  much  greater  loss  of  NTO  in  the  vadose  zone  is  due  to  the  longer 
transit  time  in  the  vadose  zone  coupled  with  a  relatively  low  HL  used  for 
NTO. 

Aquifer  concentrations  of  NTO  for  this  site  are  weakly  sensitive  to  NTO  HL 
in  surface  soil  and  mildly  sensitive  to  NTO  HL  in  aquifer,  but  are  strongly 
affected  by  HL  in  the  vadose  zone.  The  high  sensitivity  to  vadose  zone  HL 
is  due  to  the  relatively  long  vadose  zone  transit  time  coupled  with  the 
range  of  recommended  HL  values. 

3.5.3  Soil  partitioning 

The  sensitivity  of  aquifer  concentrations  to  NTO  Kd  values  in  each  medium 
was  evaluated.  This  sensitivity  analysis  was  conducted  by  varying  Kd  in  one 
medium  at  a  time  while  setting  the  Kd  values  to  the  recommended  means 
for  the  other  two  media  and  then  comparing  computed  aquifer  peak  well 
concentration  to  the  peak  computed  with  recommended  mean  Kd  values 
for  all  three  media  (see  Figure  2).  All  HL  values  were  set  to  too  years  for 
these  model  runs. 

Using  the  recommended  maximum  NTO  Kd  value  in  surface  soil  of  0.91 
L/kg  (see  Table  3)  resulted  in  an  aquifer  peak  well  concentration  of  4.17 
ppb  compared  to  4.29  ppb  for  the  recommended  mean  NTO  Kd  value  in 
surface  soil  of  0.52  L/kg,  which  is  a  reduction  of  2.8  %.  Using  the  recom¬ 
mended  minimum  NTO  Kd  value  in  surface  soil  of  0.36  L/kg  (see  Table  3) 
resulted  in  an  aquifer  peak  well  concentration  of  4.47  ppb  compared  to 
4.29  ppb  for  the  recommended  mean  NTO  Kd  value,  or  an  increase  of  4.2 
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%.  Thus,  aquifer  concentrations  of  NTO  for  this  site  have  a  low  sensitivity 
to  surface  soil  NTO  Kd  values.  However,  it  is  recognized  that  the  sensitivity 
to  surface  soil  Kd  values  would  increase  for  lower  HL  in  soil. 

The  maximum  and  minimum  NTO  Kd  values  were  applied  for  the  aquifer 
with  all  other  Kd  values  for  the  other  two  media  set  to  recommended  mean 
values  (see  Table  3).  These  simulations  showed  that  aquifer  concentrations 
of  NTO  were  quite  close  (less  than  0.01  ppb  difference)  to  those  computed 
using  the  recommended  mean  Kd  value.  Thus,  aquifer  concentrations  of 
NTO  for  this  site  are  practically  insensitive  to  aquifer  NTO  Kd  values. 

Using  the  maximum  recommended  value  for  NTO  vadose  Kd  of  0.64  L/kg, 
the  peak  well  concentration  for  NTO  was  2.51  ppb  (41  %  decrease  from 
4.29  ppb)  and  occurred  after  68  years  rather  than  47  years  using  Kd  of 
0.36  L/kg.  Using  the  minimum  recommended  Kd  value  of  0.25  for  the  va¬ 
dose  zone,  a  peak  well  concentration  for  NTO  of  5.44  ppb  (27  %  increase 
from  4.29  ppb)  was  computed  after  38  years.  The  percent  of  changes  are 
disproportional  due  to  the  increased  effect  of  degradation  for  the  longer 
transit  time  associated  with  a  higher  Kd  value.  The  HL  for  all  media  was  set 
to  too  years.  It  is  concluded  that  aquifer  NTO  concentrations  are  moder¬ 
ately  sensitive  to  Kd  values  in  the  vadose  zone  for  this  site.  The  sensitivity 
of  NTO  Kd  value  in  the  vadose  zone  depends  on  NTO  HL  in  the  vadose 
zone,  where  lower  vadose  HL  increases  sensitivity  of  vadose  Kd  value. 

3.5.4  Conclusions 

Aquifer  NTO  concentrations  for  this  site  are  practically  insensitive  to  Kd 
values  in  soil  and  aquifer,  but  are  moderately  sensitive  to  Kd  values  in  the 
vadose  zone.  However,  assuming  a  very  long  HL,  aquifer  concentrations 
vary  with  vadose  zone  Kd  due  to  the  fact  that  the  soil  NTO  loading  imposed 
at  this  site  extended  for  a  relatively  brief  finite  period  of  time  (10  years). 
Had  the  loadings  been  steady  state,  steady-state  aquifer  concentrations 
would  eventually  be  reached,  and  for  a  very  long  HL,  Kd  values  would  have 
had  little  effect  on  the  steady-state  concentration  but  would  have  affected 
the  time  required  to  reach  steady  state.  Aquifer  concentrations  of  NTO  for 
this  site  are  weakly  sensitive  to  NTO  HL  in  surface  soil  and  mildly  sensi¬ 
tive  to  NTO  HL  in  aquifer,  but  are  strongly  affected  by  HL  in  the  vadose 


zone. 
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3.6  NQ  sensitivity 

3.6.1  Residence  and  transit  times 

A  simulation  of  NQ  was  conducted  to  evaluate  its  residence  time  in  surface 
soil  and  transit  times  through  the  vadose  zone  and  aquifer  from  vadose 
entry  to  receptor  well.  This  simulation  was  conducted  in  the  same  manner 
as  described  above  for  a  general  tracer,  DNAN,  and  NTO.  Thus,  NQ  was 
treated  as  a  miscible  constituent  with  a  constant  soil  loading  rate  of  50 
years  and  no  degradation  in  any  medium.  The  exception  was  that  the  Kd 
values  were  set  to  the  recommended  mean  values  for  NQ,  or  0.3  L/kg  for 
all  media  (Table  3). 

The  results  of  this  simulation  showed  that  the  residence  time  in  soil  is 
about  2  months,  and  the  transit  time  through  the  vadose  zone  and  from 
vadose  to  the  receiving  well  are  about  35  years  and  4  months,  respectively. 
NQ  moves  through  each  medium  faster  than  DNAN  and  NTO  due  to  the 
lower  Kd  values  of  NQ. 

3.6.2  Half-lives 

The  sensitivity  of  aquifer  concentrations  to  half-life  of  NQ  in  each  medium 
was  evaluated.  This  sensitivity  analysis  was  conducted  by  varying  the  HL 
in  one  medium  at  a  time  while  setting  the  HL  to  too  years  for  the  other 
two  media  and  then  comparing  computed  aquifer  peak  well  concentration 
to  the  peak  computed  with  HL  set  to  too  years  for  all  three  media  (see 
Figure  2).  All  Kd  values  were  set  to  baseline  values  for  these  model  runs 
(see  recommended  mean  values  in  Table  3). 

Using  the  maximum  recommended  NQ  HL  of  3  years  in  soil  (see  Table  5) 
and  too  years  below  ground  in  vadose  and  aquifer  only  slightly  affected 
the  aquifer  well  concentration,  reducing  the  peak  concentration  from  5.05 
ppb  for  100-year  half-life  in  all  media  to  4.75  ppb,  or  5.9  %  reduction.  Us¬ 
ing  the  minimum  recommended  NQ  half-life  in  soil  of  0.3  years  and  100 
years  below  ground  resulted  in  reducing  well  peak  concentration  to  4.14 
ppb,  or  18  %  reduction  from  results  with  a  100-year  half-life  in  all  media. 
Thus,  aquifer  concentration  of  NQ  is  mildly  sensitive  to  NQ  HL  in  soil  for 
the  range  of  recommended  NQ  HL  values. 

With  the  maximum  recommended  NQ  HL  of  3  years  in  aquifer  (see  Table 
5)  and  100  years  in  soil  and  vadose  zone,  the  aquifer  peak  well  concentra- 
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tion  was  reduced  from  5.05  ppb  for  100-year  half-life  in  all  media  to  4.58 
ppb,  or  a  reduction  of  9.3  %.  Using  the  minimum  recommended  half-life 
in  aquifer  of  0.3  year  and  100  years  for  soil  and  vadose  zone  resulted  in  a 
peak  well  concentration  for  NQ  of  1.9  ppb,  or  a  62  %  reduction  compared 
with  using  an  HL  of  100  years  in  all  three  media.  Aquifer  NQ  concentra¬ 
tions  are  moderately  sensitive  to  the  HL  of  NQ  in  the  aquifer  for  the  range 
of  recommended  NQ  HL  values. 

Using  a  NQ  HL  of  3  years  in  the  vadose  zone  and  100  years  in  the  other 
two  media  resulted  in  essentially  zero  aquifer  concentration  of  NQ.  The 
much  greater  loss  of  NQ  in  the  vadose  zone  is  due  to  the  longer  transit 
time  in  the  vadose  zone  coupled  with  a  relatively  low  HL  for  NQ. 

Aquifer  concentrations  of  NQ  for  this  site  are  weakly  sensitive  to  NQ  HL  in 
surface  soil  and  moderately  sensitive  to  NQ  HL  in  aquifer,  but  are  strongly 
affected  by  HL  in  the  vadose  zone.  The  high  sensitivity  to  vadose  zone  HL 
is  due  to  the  relatively  long  vadose  zone  transit  time  coupled  with  the 
range  of  recommended  HL  values. 

3.6.3  Soil  partitioning 

The  sensitivity  of  aquifer  concentrations  to  NQ  Kd  values  in  each  medium 
was  evaluated.  This  sensitivity  analysis  was  conducted  by  varying  Kd  in  one 
medium  at  a  time  while  setting  the  Kd  values  to  the  recommended  means 
for  the  other  two  media  and  then  comparing  computed  aquifer  peak  well 
concentration  to  the  peak  computed  with  recommended  mean  Kd  values 
for  all  three  media  (see  Figure  2).  All  HL  values  were  set  to  too  years  for 
these  model  runs. 

Using  the  recommended  maximum  and  minimum  NQ  Kd  value  in  surface 
soil  of  0.6  and  0.1  L/kg  (see  Table  3)  resulted  in  aquifer  peak  well  concen¬ 
trations  very  close  (less  than  0.02  ppb)  to  the  value  of  5.05  ppb,  which  was 
computed  using  the  recommended  mean  NQ  Kd  value  in  surface  soil  of  0.3 
L/kg.  Thus,  aquifer  concentrations  of  NQ  for  this  site  are  practically  insen¬ 
sitive  to  surface  soil  NQ  Kd  values.  However,  it  is  recognized  that  the  sensi¬ 
tivity  to  surface  soil  Kd  values  would  increase  for  lower  HL  in  soil,  where  a 
value  of  too  years  was  used  for  this  analysis. 

The  maximum  and  minimum  NQ  Kd  values  were  applied  for  the  aquifer 
with  all  other  Kd  values  for  the  other  two  media  set  to  recommended  mean 
values  (see  Table  3).  These  simulations  showed  that  aquifer  concentrations 
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of  NQ  were  quite  close  (less  than  0.03  ppb  difference)  to  those  computed 
using  the  recommended  mean  Kd  value.  Thus,  aquifer  concentrations  of 
NQ  for  this  site  are  practically  insensitive  to  aquifer  NQ  Kd  values. 

Using  the  maximum  recommended  value  for  NQ  vadose  Kd  of  0.6  L/kg, 
the  peak  well  concentration  for  NQ  was  2.8  ppb  (45  %  decrease  from  5.05 
ppb)  and  occurred  after  65  years  rather  than  41  years  using  Kd  of  0.3  L/kg. 
Using  the  minimum  recommended  Kd  value  of  0.1  for  the  vadose  zone,  a 
peak  well  concentration  for  NQ  of  7.49  ppb  (48  %  increase  from  5.05  ppb) 
was  computed  after  24  years.  It  is  concluded  that  aquifer  NQ  concentra¬ 
tions  are  moderately  sensitive  to  Kd  values  in  the  vadose  zone  for  this  site. 
The  sensitivity  of  NQ  Kd  value  in  the  vadose  zone  depends  on  NQ  HL  in 
the  vadose  zone.  However,  the  sensitivity  of  NQ  vadose  Kd  value  is  of  little 
consequence  for  the  NQ  HL  values  shown  in  Table  5  that  vary  between  .3 
and  3  years  since  an  HL  of  3  years  resulted  in  low  NQ  aquifer  concentra¬ 
tions  (0.15  ppb  with  an  NQ  vadose  Kd  value  of  0.1  L/kg). 

3.6.4  Conclusions 

Aquifer  NQ  concentrations  for  this  site  are  practically  insensitive  to  Kd 
values  in  soil  and  aquifer,  but  are  moderately  sensitive  to  Kd  values  in  the 
vadose  zone.  However,  assuming  a  very  long  HL,  aquifer  concentrations 
vary  with  vadose  zone  Kd  due  to  the  fact  that  the  soil  NQ  loading  imposed 
at  this  site  extended  for  a  relatively  brief  finite  period  of  time  (10  years). 
Had  the  loadings  been  steady  state,  steady-state  aquifer  concentrations 
would  eventually  be  reached,  and  for  a  very  long  HL,  Kd  values  would  have 
little  effect  on  the  steady-state  concentration  but  would  affect  the  time  to 
reach  steady  state.  Aquifer  concentrations  of  NQ  for  this  site  are  weakly 
sensitive  to  NQ  HL  in  surface  soil  and  moderately  sensitive  to  NQ  HL  in 
aquifer,  but  are  strongly  affected  by  HL  in  the  vadose  zone.  The  high  sensi¬ 
tivity  to  vadose  zone  HL  is  due  to  the  relatively  long  vadose  zone  transit 
time  coupled  with  the  range  of  recommended  HL  values. 

3.7  Summary 

It  is  important  to  know  soil  residence  time  and  transit  time  through  the 
vadose  zone  and  aquifer  to  more  fully  understand  why  some  model  inputs 
are  more  important  than  others.  The  amount  of  constituent  degradation  in 
various  media  depends  on  the  degradation  rate  (i.e.,  HL)  in  a  medium  and 
the  time  spent  in  that  medium.  The  time  spent  in  a  medium  is  affected  by 
the  soil  partitioning  coefficient.  The  tracer  simulation  for  this  study  site 
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showed  that  the  soil  residence  time  and  transit  time  through  the  aquifer 
(from  vadose  entry  to  receptor  well)  is  only  a  few  months,  whereas  the 
transit  time  through  the  vadose  zone  is  approximately  14  years  for  a  non¬ 
adsorbing  tracer.  With  soil  partitioning  of  the  EC,  the  residence  and  trans¬ 
it  times  are  extended;  i.e.,  transport  is  retarded.  The  vadose  zone  is  by  far 
the  most  important  medium  for  this  site  with  regard  to  EC  fate. 

The  sensitivity  results  obtained  through  the  model  applications  for  the 
three  EC  to  Demo  Area  2,  MMR,  are  summarized  in  Table  6.  It  should  be 
recognized  that  the  high  sensitivity  of  some  inputs  are  related  to  the  site 
characteristics.  For  example,  this  site  has  a  much  longer  transit  time 
through  the  vadose  than  through  the  aquifer,  but  the  reverse  can  be  the 
case  for  other  sites. 


Table  6.  Sensitivity  of  aquifer  EC  concentrations  to  inputs  for  soil  partitioning  and 

degradation. 


EC 

Soil  HL 

Aquifer  HL 

Vadose  HL 

Soil  Ka 

Aquifer  Kd 

Vadose  Kd 

DNAN 

Moderate 

High 

Very  high 

Low 

Very  low 

Moderate 
to  high 

NTO 

Low 

Low  to 
moderate 

High  to 
very  high 

Low 

Very  low 

Moderate 

NQ 

Low  to 
moderate 

Moderate 

Very  High 

Very  low 

Very  low 

Moderate 

Examination  of  Table  6  reveals  that  model  results  are  either  highly  or  very 
highly  sensitive  to  model  input  values  for  vadose  zone  HL.  Model  results 
range  from  low  to  moderate  sensitivity  to  soil  HL  input  values.  Model  re¬ 
sults  range  from  low  to  high  sensitivity  to  aquifer  HL  input  values.  Model 
results  are  either  very  lowly  or  lowly  sensitive  to  model  input  values  for 
soil  and  aquifer  Kd  values.  Results  are  moderately  to  highly  sensitive  to  va¬ 
dose  Kd  values.  Therefore,  it  is  concluded  that  estimation  of  HL,  particu¬ 
larly  for  vadose  zone  and  groundwater,  is  critically  important  for  the 
accuracy  of  model  predictions.  Furthermore,  accurate  estimation  of  HL  is 
considerably  more  important  than  accurate  estimation  of  Kd.  However,  it 
is  recognized  that  Kd  can  affect  the  influence  of  HL  in  each  medium,  where 
high  Kd  values  cause  HL  to  exert  more  influence  on  model  results. 
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4  Artillery  Impact  Area  (AIA),  USMA, 
Application 

4.1  Background  for  model  inputs 

This  chapter  presents  the  results  of  application  of  TREECS™  to  the  AIA  at 
USMA  to  assess  the  fate  of  the  three  EC  in  surface  water  as  affected  by  un¬ 
certain  inputs  for  soil-water  partitioning  and  half-life.  The  site  description 
and  model  inputs  for  this  study  site  were  described  previously  (Dortch 
2012  for  overall  USMA  ranges  and  Dortch  et  al.  2017  for  AIA  application) 
so  are  not  repeated  here.  However,  basic  information  and  key  features  of 
the  model  inputs  regarding  this  application  are  described  below  to  provide 
the  background  required  for  discussing  model  results. 

Two  modeling  components  of  TREECS™  were  used,  which  included  the 
Tier  2  soil  model  and  the  Contaminant  Model  for  Streams  (CMS).  The  op¬ 
tion  for  average  annual  hydrology,  rather  than  daily  varying  hydrology, 
was  used  in  the  soil  model  for  this  application  to  facilitate  long-term  (cen¬ 
tury)  application.  A  constant  EC  residue  loading  rate  of  49.8  kilograms  per 
year  was  applied  for  100  years.  This  loading  rate  is  the  same  as  the  HE 
RDX  loading  rate  used  in  the  previous  applications  to  this  site  (Dortch  et 
al.  2017);  the  rate  was  calculated  with  the  TREECS™  impact  zone  MC  res¬ 
idue  calculator  based  on  the  numbers  and  types  of  munitions  fired  on 
these  artillery  ranges  as  described  by  Dortch  et  al.  (2017). 

The  model  was  run  for  100  years  using  the  adaptive  time  step  feature  and 
minimum  time  step  of  0.001  year.  Zero  initial  soil  contamination  of  EC 
was  assumed.  All  of  the  soil  infiltrating  water  and  EC  mass  flux  was  as¬ 
sumed  to  be  returned  to  surface  flow  via  soil  interflow  due  to  bedrock  that 
underlies  surface  soil,  thus  preventing  percolation.  Solid  and  adsorbed 
phase  EC  were  assumed  to  be  nondegradable;  thus,  a  high  value  of  1E20 
years  was  used  in  the  soil  model  for  half-life  of  solid-phase  EC  and  EC  ad¬ 
sorbed  to  soil,  and  a  degradation  rate  of  zero  was  used  in  CMS  forEC  ad¬ 
sorbed  to  stream  benthic  and  suspended  sediments.  Soil  properties  for 
sandy  loam  were  used  for  range  impact  zone  surface  soil.  Chemical- 
specific  property  inputs  of  the  EC  are  presented  in  Table  1  and  Table  2. 

The  mean  and  range  of  recommended  EC  soil  and  sediment  Kd  values  are 
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presented  in  Table  3  and  Table  4,  respectively.  The  mean  and  range  of  rec¬ 
ommended  EC  half-lives  are  presented  in  Table  5. 

The  initial  particle  size  of  solid  phase  EC  residue  was  set  to  5,401  pm  using 
an  approach  recommended  by  Dortch  et  al.  (2017),  which  was  based  on 
modeling  of  IMX-101  fate  laboratory  studies.  This  approach  is  described  in 
the  previous  chapter.  The  only  difference  in  the  calculation  of  particle  size 
here  from  the  previous  calculation  for  Demo  Area  2  is  that  the  representa¬ 
tive  size  of  RDX  was  12,000  pm  rather  than  6,000  used  for  Demo  Area  2. 
The  RDX  particle  size  was  assumed  to  correspond  to  the  DNAN  particle 
size  with  the  procedure.  The  particle  diameter  of  5,401  pm  was  used  as  the 
input  value  for  each  of  the  three  EC. 

The  target  receptor  location  is  Popolopen  Brook,  which  receives  drainage 
from  the  area  of  interest  (AOI),  or  the  AIA  for  the  USMA  artillery  ranges. 
For  the  previous  AIA  application  (Dortch  et  al.  2017),  the  target  receptor 
location  was  set  at  50  m  downstream  of  the  AIA  influent  into  Popolopen 
Brook  (the  beginning  of  the  modeled  stream  reach),  which  is  the  approxi¬ 
mate  location  where  stream  samples  were  collected.  Thus,  this  location 
was  also  used  for  the  present  application.  However,  a  second  target  recep¬ 
tor  location  further  downstream  was  also  used  to  evaluate  the  importance 
of  EC  degradation  within  the  stream.  The  second  location  was  arbitrarily 
placed  5,000  m  downstream  from  the  beginning  of  the  modeled  stream 
reach.  This  much  longer  distance  was  chosen  to  provide  ample  travel  time 
for  degradation  to  have  a  perceivable  effect.  All  other  inputs  pertaining  to 
meteorology,  soil  conditions,  hydrology,  stream  conditions,  and  hydraulics 
are  provided  by  Dortch  et  al.  (2017). 

4.2  Baseline  results 

TREECS™  was  applied  using  the  inputs  referenced  and  discussed  in  the 
previous  section  to  establish  baseline  results  for  each  EC.  The  recom¬ 
mended  mean  soil  Kd  values  and  mean  HL  for  surface  soil,  stream  sedi¬ 
ments,  and  surface  water  were  used  for  the  baseline  application.  RDX  was 
also  included  in  the  baseline  results  for  comparing  the  EC.  The  inputs  for 
RDX  were  the  same  as  those  of  the  validation  application  presented  by 
Dortch  et  al.  (2017).  The  HL  of  RDX  in  AOI  soil  was  set  to  10  years,  and 
zero  degradation  rate  of  RDX  was  used  in  the  stream  for  the  validation  ap¬ 
plication. 
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The  stream  total  concentrations  (particulate  plus  dissolved)  in  parts  per 
billion  (ppb),  or  micrograms  per  liter,  for  the  three  EC  and  RDX  are  plot¬ 
ted  versus  time  for  the  upstream  receptor  target  location  as  shown  in  Fig¬ 
ure  3.  The  concentrations  of  the  three  EC  are  greater  than  that  of  RDX  and 
reach  a  steady-state  concentration  much  faster.  Eventually  RDX  would  al¬ 
so  reach  a  steady-state  concentration  of  about  the  same  concentration  as 
the  EC  if  the  model  was  run  for  a  longer  period  of  time.  The  much  lower 
solubility  and  the  larger  initial  particle  size  (12,000  pm)  of  RDX  are  the 
reasons  that  it  takes  so  much  longer  for  RDX  to  reach  steady  state.  Steady- 
state  concentrations  are  the  result  of  the  constant  AOI  loading  rate  of  EC 
and  RDX.  The  differences  in  steady-state  EC  concentrations  are  caused  by 
the  differences  in  soil  degradation  rates  or  HL,  where  higher  steady-state 
concentration  corresponds  with  higher  HL  in  soil. 

Figure  3.  Computed  stream  concentration  of  EC  compared  with  RDX  50  m 
downstream  of  the  AIA,  USMA,  for  baseline  conditions. 
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The  application  with  baseline  conditions  was  repeated,  except  that  the  re¬ 
ceptor  location  was  moved  further  downstream  to  5,000  m  downstream  of 
the  AIA.  The  results  of  this  application  are  shown  in  Figure  4.  With  an  HL 
of  only  one  day  for  each  EC,  a  substantial  reduction  in  steady-state  con¬ 
centrations  (compared  with  results  in  Figure  3)  occurs  over  the  longer 
stream  reach  where  a  longer  travel  time  allows  more  degradation.  The 
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concentrations  of  RDX  are  essentially  the  same  in  Figure  3  and  Figure  4 
due  to  the  assumption  of  no  RDX  degradation  within  the  stream. 


Figure  4.  Computed  stream  concentration  of  EC  compared  with  RDX  5000  m 
downstream  of  the  AIA,  USMA,  for  baseline  conditions. 


It  is  noted  that  the  above  baseline  results  are  not  totally  realistic  for  actual 
field  conditions.  Stream  concentrations  of  MC  in  the  field  have  been  found 
to  be  highly  transient,  varying  from  below  detection  to  levels  of  several 
parts  per  billion.  The  higher  concentrations  occur  during  and  immediately 
following  significant  rainfall  events,  and  return  to  below  detection  as  rain¬ 
fall  and  runoff  cease  (Dortch  et  al.  2017).  Thus,  continuous,  steady-state 
concentrations  as  shown  above  are  not  expected  to  occur  in  the  field,  but 
they  serve  the  purposes  of  this  study. 

Tracer  modeling  was  not  performed  for  the  AIA  study  area  since  travel 
time  within  the  stream  is  relatively  short  and  can  be  easily  estimated. 

Short  travel  times  in  surface  water  result  in  less  importance  for  accurately 
estimating  surface  water  degradation  rates  or  HL.  The  sensitivity  of  uncer¬ 
tain  Kd  and  HL  values  for  this  surface  water  system  are  discussed  in  the 
sections  below. 
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4.3  DNAN  sensitivity 

4.3.1  Half-lives 

Given  the  highly  transient  nature  of  expected  stream  EC  concentrations, 
the  short  travel  times  in  streams,  and  the  expected  high  photodegradation 
rates  yielding  HL  on  the  order  of  one  day,  HL  sensitivity  was  focused  pre¬ 
dominantly  on  soil  HL  values.  All  Kd  values  were  set  to  baseline  values  for 
these  model  runs  (see  recommended  mean  values  in  Table  3).  All  other  in¬ 
puts  were  the  same  as  the  baseline  conditions  with  the  exception  of  the 
length  of  simulation,  which  was  shortened  to  50  years  since  steady-state 
conditions  are  reached  relatively  quickly. 

Monte  Carlo  uncertainty  analysis  was  conducted  by  varying  only  the  soil 
HL  values  and  monitoring  stream  total  EC  concentrations  at  the  upstream 
station  (50  m  downstream  of  the  AIA).  The  mean  soil  HL  for  DNAN  is  0.5 
years  (see  Table  5)  with  lower  and  upper  bounds  of  0.02  and  2.2  years.  A 
log  normal  distribution  was  assumed  for  the  uncertainty  analysis.  The  HL 
values  for  DNAN  in  Table  5  were  converted  to  degradation  rates  as  re¬ 
quired  by  the  Monte  Carlo  Sensitivity/Uncertainty  (S/U)  module  by  as¬ 
suming  a  first  order  decay.  Thus,  the  upper  and  lower  bounds  on 
degradation  rate  were  set  to  34.7  and  0.32  per  year,  respectively,  within 
the  S/U  module  UI.  The  log  of  the  mean  degradation  rate  was  set  to  0.14 
per  year  using  the  HL  of  0.5  year  and  was  entered  into  the  S/U  UI.  The  dif¬ 
ference  in  the  logs  of  the  upper  and  lower  degradation  rates  was  divided  by 
6  to  estimate  the  standard  deviation  log  of  0.34  per  year  that  was  input  in 
the  S/U  UI.  The  S/U  module  was  executed  for  50  iterations,  which  was 
enough  to  reach  consistent  results  with  the  Latin  Hypercube  sampling 
procedure. 

The  results  of  the  above  S/U  application  are  plotted  in  Figure  5  for  base¬ 
line  results  and  upper  and  lower  95  %  confidence  bands  associated  with 
uncertainty  of  soil  DNAN  HL.  This  figure  shows  that  stream  DNAN  con¬ 
centrations  are  substantially  affected  by  soil  HL.  Two  deterministic  runs 
were  executed  with  the  soil  HL  set  to  the  upper  and  lower  bounds  of  2.2 
and  0.02  years.  The  resulting  steady-state  stream  concentrations  of  DNAN 
were  3.08  and  0.57  ppb,  respectively.  The  concentration  of  3.08  ppb  is 
slightly  greater  than  the  upper  confidence  limit  (UCL),  and  the  concentra¬ 
tion  of  0.57  ppb  is  about  a  fourth  of  the  lower  confidence  limit  (LCL).  This 
rather  large  skewness  of  the  confidence  bands  is  due  to  using  a  log  normal 
distribution  with  a  mean  HL  of  0.5  years. 
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Figure  5.  Computed  stream  concentration  of  DNAN  50  m  downstream  of  the  AIA, 
USMA,  with  95  %  confidence  bounds  for  uncertainty  of 
soil  HL  compared  to  baseline  conditions. 
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An  application  was  performed  to  examine  the  sensitivity  of  the  DNAN 
stream  degradation  rate  for  the  second  receptor  location  of  5,000  m 
downstream  of  the  AIA.  The  stream  degradation  rate  was  set  to  zero  for 
DNAN  in  this  application.  The  resulting  steady-state  stream  concentration 
of  DNAN  at  the  receptor  location  for  this  application  was  2.71  ppb,  which 
is  about  the  same  as  the  steady-state  DNAN  concentration  at  the  upstream 
receptor  location  using  a  degradation  rate  of  1.0  per  day  (Figure  3).  Figure 
4  shows  that  the  steady-state  DNAN  concentration  at  the  downstream  re¬ 
ceptor  location  with  a  degradation  rate  of  1.0  per  day  is  2.12  ppb.  Thus,  the 
rather  high  degradation  rate  of  1.0  per  day  results  in  only  a  0.59  ppb  de¬ 
crease  in  stream  concentration  over  5,000  m  compared  with  no  degrada¬ 
tion.  It  is  concluded  that  stream  concentrations  of  DNAN  are  more 
sensitive  to  soil  HL  than  stream  degradation  rate.  This  higher  sensitivity  is 
due  to  the  longer  DNAN  residence  time  in  soil  compared  to  travel  time  in 
the  stream. 

4.3.2  Soil  partitioning 

The  sensitivity  of  stream  DNAN  concentrations  to  DNAN  Kd  values  in  soil 
was  evaluated.  Monte  Carlo  uncertainty  analysis  was  conducted  by  varying 
only  the  soil  DNAN  Kd  values  and  monitoring  stream  total  DNAN  concen- 
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trations  at  the  upstream  station  (50  m  downstream  of  the  AIA).  A  normal 
distribution  was  assumed  with  a  mean  soil  DNAN  Kd  value  of  2.98  L/kg 
and  lower  and  upper  bounds  of  1.33  and  4.66  L/kg  (see  Table  3).  The  dif¬ 
ference  in  the  bound  values  divided  by  6  resulted  in  a  standard  deviation 
of  0.56  L/kg.  All  four  values  were  entered  into  the  S/U  UI,  and  the  S/U 
module  was  run  for  50  iterations. 


The  results  of  the  above  S/U  application  are  plotted  in  Figure  6  for  base¬ 
line  results  and  upper  and  lower  95  %  confidence  band  associated  with  un¬ 
certainty  of  soil  DNAN  Kd  values.  This  figure  shows  that  stream  DNAN 
concentrations  are  hardly  affected  by  soil  Kd  values.  Two  deterministic 
runs  were  executed  with  the  soil  Kd  values  set  to  the  upper  and  lower 
bounds  of  4.66  and  1.33  L/kg.  The  resulting  steady-state  stream  concen¬ 
trations  of  DNAN  were  2.72  and  2.70  ppb,  respectively,  confirming  that 
soil  DNAN  Kd  values  have  very  little  effect  on  stream  concentrations  of 
DNAN. 


Figure  6.  Computed  stream  concentration  of  DNAN  50  m  downstream  of  the  AIA, 
USMA,  with  95  %  confidence  bounds  for  uncertainty  of 
soil  Kd  values  compared  to  baseline  conditions. 


Surface  water  total  concentration  in  Popolopen  Brook 
downstream  of  AIA,  USMA 


The  sensitivity  of  DNAN  Kd  values  in  water  and  sediments  was  evaluated 
by  applying  the  baseline  conditions  with  the  maximum  and  minimum  val¬ 
ues  for  DNAN  shown  in  Table  4.  The  same  sediment  Kd  value  was  used  for 
both  benthic  and  suspended  sediments.  The  computed  steady-state  stream 
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DNAN  concentrations  were  the  same  for  both  Kd  values,  indicating  insen¬ 
sitivity  to  varying  stream  DNAN  Kd  values. 

4.3.3  Conclusions 

Surface  water  stream  DNAN  concentrations  are  moderately  sensitive  to 
soil  HL  and  less  sensitive  to  stream  degradation  rate.  Stream  DNAN  con¬ 
centrations  for  this  site  are  insensitive  to  Kd  values  in  soil  and  stream. 

4.4  NTO  sensitivity 

4.4.1  Half-lives 

Monte  Carlo  uncertainty  analysis  was  conducted  by  varying  only  the  soil 
NTO  HL  values  and  monitoring  stream  total  NTO  concentrations  at  the 
upstream  station  (50  m  downstream  of  the  AIA).  The  mean  soil  HL  for 
NTO  is  3.0  years  (see  Table  5)  with  lower  and  upper  bounds  of  1.0  and 
10.0  years.  A  log  normal  distribution  was  assumed  for  the  uncertainty 
analysis.  The  HL  values  for  NTO  in  Table  5  were  converted  to  degradation 
rates  as  required  by  the  Monte  Carlo  S/U  module  by  assuming  a  first  order 
decay.  Thus,  the  upper  and  lower  bounds  on  the  degradation  rate  were  set 
to  0.69  and  0.07  per  year,  respectively,  within  the  S/U  module  UI.  The  log 
of  the  mean  degradation  rate  was  set  to  -0.636  per  year  using  the  HL  of 
3.0  years  and  was  entered  into  the  S/U  UI.  The  difference  in  the  logs  of  the 
upper  and  lower  degradation  rates  was  divided  by  six  to  estimate  the 
standard  deviation  log  of  0.167  per  year  that  was  input  in  the  S/U  UI.  The 
S/U  module  was  executed  for  50  iterations,  which  was  enough  to  reach 
consistent  results  with  the  Latin  Hypercube  sampling  procedure. 

The  results  of  the  above  S/U  application  are  plotted  in  Figure  7  for  base¬ 
line  results  and  the  upper  and  lower  95  %  confidence  band  associated  with 
uncertainty  of  soil  NTO  HL.  This  figure  shows  that  stream  NTO  concentra¬ 
tions  are  hardly  affected  by  soil  HL.  Two  deterministic  runs  were  executed 
with  the  soil  HL  set  to  the  upper  and  lower  bounds  of  10  and  1  years.  The 
resulting  steady-state  stream  concentrations  of  NTO  were  3.19  and  2.93 
ppb,  respectively,  which  confirmed  the  low  sensitivity  of  stream  NTO  con¬ 
centrations  to  soil  NTO  HL.  Another  S/U  run  was  made  assuming  a  nor¬ 
mal  distribution  with  the  mean  and  bounds  set  the  same  as  described 
above.  The  results  of  this  run  were  nearly  identical  to  those  shown  in  Fig¬ 
ure  7.  NTO  has  a  lower  residence  time  in  soil  than  DNAN  due  to  its  lower 
soil  partitioning,  thus  rendering  soil  HL  as  less  sensitive. 
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Figure  7.  Computed  stream  concentration  of  NTO  50  m  downstream 
of  the  AIA,  USMA,  with  95  %  confidence  bounds  for  uncertainty 
of  soil  HL  compared  to  baseline  conditions. 
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The  sensitivity  of  stream  NTO  degradation  rate  at  the  second  receptor  lo¬ 
cation  of  5,000  m  downstream  of  the  AIA  is  very  similar  to  that  of  DNAN. 
The  rather  high  degradation  rate  of  l.o  per  day  results  in  a  0.68  ppb  de¬ 
crease  in  steady-state  stream  concentration  of  NTO  over  5,000  m  com¬ 
pared  with  no  degradation.  It  is  concluded  that  stream  concentrations  of 
NTO  are  more  sensitive  to  stream  degradation  rate  than  soil  HL  of  NTO. 
The  low  sensitivity  of  stream  concentrations  to  soil  HL  is  due  to  a  shorter 
soil  residence  time  (compared  to  DNAN)  coupled  with  higher  soil  HL  of 
NTO.  The  stream  travel  time  of  NTO  is  the  same  as  DNAN.  The  very  wide 
range  of  stream  degradation  rate  from  none  to  one  per  day  is  the  major 
driver  resulting  in  significant  sensitivity  to  stream  degradation  rate. 

4.4.2  Soil  partitioning 

The  sensitivity  of  stream  NTO  concentrations  to  NTO  Kd  values  in  soil  was 
evaluated.  Monte  Carlo  uncertainty  analysis  was  conducted  by  varying  on¬ 
ly  the  soil  NTO  Kd  values  and  monitoring  stream  total  NTO  concentrations 
at  the  upstream  station  (50  m  downstream  of  the  AIA).  A  normal  distribu¬ 
tion  was  assumed  with  a  mean  soil  NTO  Kd  value  of  0.41  L/kg  and  upper 
and  lower  bounds  of  0.72  and  0.29  L/kg  (see  Table  3).  The  difference  in 
the  bound  values  divided  by  six  resulted  in  a  standard  deviation  of  0.072 
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L/kg.  All  four  values  were  entered  into  the  S/U  UI,  and  the  S/U  module 
was  run  for  50  iterations. 

The  results  of  the  above  S/U  application  look  similar  to  the  results  plotted 
in  Figure  7;  thus  exhibiting  almost  no  sensitivity  to  soil  NTO  Kd  values. 

The  sensitivity  of  Kd  values  of  NTO  in  water  and  sediments  was  evaluated 
by  applying  the  baseline  conditions  with  the  maximum  and  minimum  val¬ 
ues  for  NTO  shown  in  Table  4.  The  computed  steady-state  stream  NTO 
concentrations  were  the  same  for  both  Kd  values,  indicating  insensitivity  to 
varying  stream  NTO  Kd  values. 

4.4.3  Conclusions 

Surface  water  stream  NTO  concentrations  are  nearly  insensitive  to  soil  HL. 
Stream  NTO  concentrations  show  about  the  same  sensitivity  to  stream 
degradation  rate  as  was  found  for  DNAN.  Stream  NTO  concentrations  for 
this  site  are  insensitive  to  NTO  Kd  values  in  soil  and  stream. 

4.5  NQ  sensitivity 

4.5.1  Half-lives 

Monte  Carlo  uncertainty  analysis  was  conducted  by  varying  only  the  soil 
NQ  HL  values  and  monitoring  stream  total  NQ  concentrations  at  the  up¬ 
stream  station  (50  m  downstream  of  the  AIA).  The  mean  soil  HL  for  NQ  is 
1.0  year  (see  Table  5)  with  lower  and  upper  bounds  of  0.3  and  3.0  years.  A 
log  normal  distribution  was  assumed  for  the  uncertainty  analysis.  The  HL 
values  for  NQ  in  Table  5  were  converted  to  degradation  rates  as  required 
by  the  Monte  Carlo  S/U  module  by  assuming  a  first  order  decay.  Thus,  the 
upper  and  lower  bounds  on  degradation  rate  were  set  to  2.31  and  0.23  per 
year,  respectively,  within  the  S/U  module  UI.  The  mean  log  degradation 
rate  was  set  to  -0.1592  per  year  using  the  HL  of  1.0  year  and  was  entered 
into  the  S/U  UI.  The  difference  in  the  logs  of  the  upper  and  lower  degrada¬ 
tion  rates  was  divided  by  6  to  estimate  the  standard  deviation  log  of  0.167 
per  year  that  was  input  in  the  S/U  UI.  The  S/U  module  was  executed  for 
50  iterations,  which  was  enough  to  reach  consistent  results  with  the  Latin 
Hypercube  sampling  procedure. 

The  results  of  the  above  S/U  application  are  plotted  in  Figure  8  for  base¬ 
line  results  and  the  upper  and  lower  95  %  confidence  band  associated  with 
uncertainty  of  soil  NQ  HL.  This  figure  shows  that  stream  NQ  concentra- 
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tions  are  slightly  sensitive  to  soil  HL  of  NQ.  Two  deterministic  runs  were 
executed  with  the  soil  HL  set  to  the  upper  and  lower  bounds  of  3  and  0.3 
years.  The  resulting  steady-state  stream  concentrations  of  NQ  were  3.12 
and  2.42  ppb,  respectively,  which  confirmed  the  fairly  low  sensitivity  of 
stream  NQ  concentrations  to  soil  NQ  HL.  NQ  responds  similarly  to  NTO 
due  to  its  high  solubility  and  relatively  low  Kd  values. 


Figure  8.  Computed  stream  concentration  of  NQ  50  m  downstream  of  the  AIA, 
USMA,  with  95  %  confidence  bounds  for  uncertainty  of 
soil  HL  compared  to  baseline  conditions. 


The  sensitivity  of  stream  NQ  degradation  rate  at  the  second  receptor  loca¬ 
tion  of  5,000  m  downstream  of  the  AIA  is  very  similar  to  that  of  DNAN 
and  NTO.  The  rather  high  degradation  rate  of  1.0  per  day  results  in  a  0.64 
ppb  decrease  in  steady-state  stream  concentration  of  NQ  over  5,000  m 
compared  with  no  degradation.  It  is  concluded  that  stream  concentrations 
of  NQ  are  more  sensitive  to  stream  degradation  rate  than  soil  HL  of  NQ. 
The  low  sensitivity  of  stream  concentrations  to  soil  HL  is  due  to  a  shorter 
soil  residence  time  (compared  to  DNAN)  coupled  with  higher  soil  HL  of 
NQ.  The  very  wide  range  of  stream  degradation  rate  from  none  to  1  per 
day  is  the  major  driver  resulting  in  higher  sensitivity  of  stream  degrada¬ 
tion  rate  compared  to  soil  HL. 
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4.5.2  Soil  partitioning 

The  sensitivity  of  stream  NQ  concentrations  to  NQ  Kd  values  in  soil  was 
evaluated.  Monte  Carlo  uncertainty  analysis  was  conducted  by  varying  on¬ 
ly  the  soil  NQ  Kd  values  and  monitoring  stream  total  NQ  concentrations  at 
the  upstream  station  (50  m  downstream  of  the  AIA).  A  normal  distribution 
was  assumed  with  a  mean  soil  NQ  Kd  value  of  0.3  L/kg  (see  Table  3)  and 
upper  and  lower  bounds  of  0.6  and  0.1  L/kg.  The  difference  in  the  bound 
values  divided  by  6  resulted  in  a  standard  deviation  of  0.083  L/kg.  All  four 
values  were  entered  into  the  S/U  UI,  and  the  S/U  module  was  run  for  50 
iterations. 

The  results  of  the  above  S/U  application  showed  practically  no  sensitivity 
to  NQ  soil  Kd  values.  The  baseline  and  95%  confidence  bands  plotted  on 
top  of  each  other.  The  sensitivity  of  Kd  values  of  NQ  in  water  and  sedi¬ 
ments  was  evaluated  by  applying  the  baseline  conditions  with  the  maxi¬ 
mum  and  minimum  values  for  NQ  shown  in  Table  4.  The  computed 
steady-state  stream  NQ  concentrations  were  the  same  for  both  Kd  values, 
indicating  insensitivity  to  varying  stream  NQ  Kd  values. 

4.5.3  Conclusions 

Surface  water  stream  NQ  concentrations  are  slightly  sensitive  to  soil  HL  of 
NQ.  Stream  NQ  concentrations  show  about  the  same  sensitivity  to  stream 
degradation  rate  as  was  found  for  DNAN  and  NTO,  with  more  sensitivity 
than  soil  HL.  Stream  NQ  concentrations  for  this  site  are  insensitive  to  NQ 
i*Cd  values  in  soil  and  stream. 

4.6  Summary 

Travel  time  in  surface  water,  especially  for  flowing  streams,  is  much  less 
important  than  transit  time  in  groundwater,  due  to  surface  water  travel 
time  being  so  short  compared  with  groundwater.  Soil  residence  time  does 
affect  receiving  stream  concentrations  of  EC  when  the  EC  degrades  more 
rapidly  in  soil.  However,  soil  residence  time  for  the  three  EC  that  were 
studied  is  relatively  short,  thus  resulting  in  lower  effect  on  reducing 
stream  concentrations  compared  with  no  degradation. 

The  sensitivity  results  obtained  through  the  model  applications  for  the 
three  EC  to  AIA,  USMA,  are  summarized  in  Table  7.  Examination  of  Table 
7  reveals  that  receiving  stream  DNAN  concentrations  are  moderately  sen- 
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sitive  to  soil  HL  for  DNAN,  but  NTO  and  NQ  stream  concentrations  are 
almost  insensitive  to  soil  HL  for  NTO  and  NQ.  This  result  is  due  to  the 
longer  residence  time  of  DNAN  in  soil  compared  with  NTO  and  NQ;  this 
condition  is  caused  by  the  higher  soil  Kd  value  of  DNAN.  Longer  soil  resi¬ 
dence  time  allows  a  greater  amount  of  degradation  in  soil,  thus  lowering 
receiving  stream  concentrations.  All  three  EC  showed  about  the  same  sen¬ 
sitivity  to  stream  degradation  rate.  The  wide  range  in  degradation  rates 
(from  zero  to  l  per  day)  that  was  tested  caused  the  moderate  level  of  sensi¬ 
tivity.  All  three  EC  exhibited  essentially  no  sensitivity  to  soil  and  stream  Kd 
values. 


Table  7.  Sensitivity  of  surface  water  stream  concentrations  for 
three  EC  to  inputs  for  degradation  and  partitioning. 


EC 

Soil  HL 

Surface 

water 

degradation 

rate 

Soil  Kd 

Surface 
water  Kd 

DNAN 

Moderate 

Moderate 

Very  low 

Very  low 

NTO 

Very  low 

Moderate 

Very  low 

Very  low 

NQ 

Low 

Moderate 

Very  low 

Very  low 

It  is  concluded  that  generally  surface  water  sites  require  much  less  accura¬ 
cy  for  estimating  soil  and  surface  water  EC  degradation  rates  than  re¬ 
quired  for  vadose  zone  and  groundwater.  However,  it  is  recognized  that 
the  importance  of  estimates  for  surface  water  degradation  rates  increases 
as  the  downstream  distance  of  (and  travel  time  to)  the  receptor  location  in 
the  receiving  stream  increases,  or  as  the  residence  time  in  the  receiving 
standing  surface  water  body  (e.g.,  ponds  and  lakes)  increases.  Also,  accu¬ 
rately  estimating  soil  degradation  rates  of  DNAN  is  more  important  than 
for  NTO  and  NQ  due  to  its  longer  soil  residence  time.  The  above  results  for 
the  AIA  site  indicate  that  for  the  expected  Kd  values  of  the  three  EC  there  is 
little  need  to  accurately  refine  the  EC  soil  and  sediment  Kd  estimates. 

Some  work  is  needed  to  verify  that  the  three  EC  Kd  values  for  sediments 
are  low  and  in  the  range  of  those  shown  in  Table  4. 
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5  Summary  and  Conclusions 

TREECS™  was  applied  to  two  study  sites  to  evaluate  the  sensitivity  of  un¬ 
certain  input  properties  for  the  fate  of  the  three  EC:  DNAN,  NTO,  and  NQ, 
which  are  components  of  the  new  insensitive  munitions  explosive  formula¬ 
tion  IMX-101.  One  study  site,  Demo  Area  2  of  MMR,  featured  groundwa¬ 
ter,  while  the  other  study  site,  AIA  of  USMA,  featured  surface  water.  Two 
types  of  input  properties  were  selected  for  a  detailed  study:  degradation 
rate  or  half-life,  HL,  and  water-to-soil  and  water-to-sediment  adsorption 
partitioning  distribution  coefficients,  or  Kd  values.  These  inputs  are  re¬ 
garded  to  be  the  most  uncertain  with  a  potentially  high  level  of  effect  on 
EC  fate.  Literature  was  reviewed  to  provide  estimates  of  the  expected  HL 
(or  degradation  rates)  and  Kd  values  for  the  three  EC,  as  well  as  the  ex¬ 
pected  range  in  these  values.  These  estimates  were  used  to  conduct  simu¬ 
lations  of  the  EC  at  each  site,  and  the  EC  fate  results  (i.e.,  concentrations 
in  receiving  groundwater  and  surface  water)  were  studied  to  draw  conclu¬ 
sions  regarding  the  sensitivity  of  results  to  these  uncertain  inputs. 

The  applications  to  Demo  Area  2,  MMR,  which  is  a  groundwater  site,  re¬ 
vealed  the  importance  of  soil  residence  time  and  transit  time  through  the 
vadose  zone  and  aquifer.  The  amount  of  constituent  degradation  in  vari¬ 
ous  media  depends  on  the  degradation  rate  (i.e.,  HL)  in  a  medium  and  the 
time  spent  in  that  medium.  The  time  spent  in  a  medium  is  affected  by  the 
soil  adsorption  partitioning  coefficient.  Adsorption  causes  transport  retar¬ 
dation.  Simulation  of  a  nondegrading  and  nonadsorbing  tracer  for  this 
study  site  showed  that  the  soil  residence  time  and  transit  time  through  the 
aquifer  (from  vadose  entry  to  receptor  well)  is  only  a  few  months,  whereas 
the  transit  time  through  the  vadose  zone  is  approximately  14  years.  Thus, 
the  vadose  zone  is  by  far  the  most  important  medium  for  this  site  with  re¬ 
gard  to  EC  fate.  However,  the  transit  time  through  the  aquifer  can  be  much 
longer  than  through  the  vadose  zone  for  other  sites  or  other  receptor  loca¬ 
tions. 

The  sensitivity  results  obtained  through  the  model  applications  for  the 
three  EC  to  Demo  Area  2,  MMR,  are  summarized  in  Table  6.  It  is  conclud¬ 
ed  that  estimation  of  HL,  particularly  for  vadose  zone  and  groundwater,  is 
critically  important  for  the  accuracy  of  model  predictions  of  receiving 
groundwater  concentrations.  Furthermore,  accurate  estimation  of  HL  is 
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considerably  more  important  than  accurate  estimation  of  Kd.  However,  it 
is  recognized  that  Kd  values  can  affect  the  amount  of  influence  of  HL, 
where  high  Kd  values  cause  HL  to  exert  more  influence  on  receiving 
groundwater  concentrations.  Thus,  the  sensitivity  of  aquifer  receiving  well 
concentrations  to  HL  in  all  media  was  greater  for  DNAN  than  for  the  other 
two  EC  due  to  the  higher  Kd  values  of  DNAN. 

The  applications  to  the  AIA,  USMA,  which  is  a  surface  water  site,  revealed 
that  travel  time  in  surface  water,  especially  for  flowing  streams,  is  much 
less  important  than  transit  time  in  groundwater,  due  to  surface  water  trav¬ 
el  time  being  so  short  compared  to  groundwater.  Soil  residence  time  does 
affect  receiving  stream  concentrations  of  EC  when  the  EC  degrades  in  soil. 
However,  soil  residence  time  for  the  three  EC  that  were  studied  is  relative¬ 
ly  short  for  this  site,  thus  resulting  in  a  modest  effect  on  reducing  stream 
concentrations  compared  with  no  degradation. 

The  sensitivity  results  obtained  through  the  model  applications  for  the 
three  EC  to  AIA,  USMA,  are  summarized  in  Table  7.  It  is  concluded  that 
generally,  surface  water  sites  require  much  less  accuracy  for  estimating 
soil  and  surface  water  EC  degradation  rates  than  what  vadose  zone  and 
groundwater  require.  However,  it  is  recognized  that  the  importance  of  es¬ 
timates  for  surface  water  degradation  rates  increases  as  the  downstream 
distance  of  (and  travel  time  to)  the  receptor  location  in  the  receiving 
stream  increases,  or  as  the  residence  time  in  the  receiving  standing  surface 
water  body  (e.g.,  ponds  and  lakes)  increases.  Also,  accurately  estimating 
soil  degradation  rates  of  DNAN  is  more  important  than  for  NTO  and  NQ, 
due  to  its  longer  soil  residence  time. 

Given  the  expected  low  Kd  values  of  the  three  EC,  there  is  minimal  need  to 
accurately  refine  the  EC  Kd  estimates,  particularly  for  surface  water  recep¬ 
tors.  Some  work  is  needed  to  verily  that  the  three  EC  Kd  values  for  sedi¬ 
ments  are  low  and  in  the  range  of  those  shown  in  Table  4. 

The  results  of  the  study  lead  to  the  conclusion  that  future  research  should 
focus  predominantly  on  determining  natural  degradation  rates  of  EC  in 
surface  soil,  vadose  zone  (partially  water-saturated  groundwater),  and  wa¬ 
ter-saturated  groundwater.  Vadose  zone  and  groundwater  should  receive 
higher  priority  surface  soil.  It  is  expected  that  EC  degradation  rates  are  not 
only  EC-specific,  but  are  also  affected  by  site-specific  media  characteris¬ 
tics,  thus  complicating  their  estimation. 
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